AD-A237  555 


NASA  Contractor  Report  187569 


ICASE 


S 


DTIC 

ELECT E 

JUN2  7  1991 

c 


SEMIANNUAL  REPORT 


October  1,  1990  through  March  31,  1991 


Contract  No.  NAS  1-18605 
May  1991 


Institute  for  Computer  Applications  in  Science  and  Engineering 
NASA  Langley  Research  Center 
Hampton,  Virginia  23665-5225 

Operated  by  the  Universities  Space  Research  Association 

BrsraiByti’oN'  stat! enhWa 

Approved  for  public  release; 

Distribution  Unlimited 

NASA 

National  Aeronautics  and 
Space  Administration 

Langley  Research  Center 

Hampton,  Virginia  23665-5225 


^'036li 

I  Him  Hill  Hill  ll/ll  mu  ■ 


02? 


CONTENTS 


Page 


Introduction . ii 

Research  in  Progress . 1 

Reports  and  Abstracts . 26 

ICASE  Interim  Reports . 47 

ICASE  Colloquia . 48 

Other  Activities . 52 


ICASE  Staff 


53 


Ac imiIm  I*r 

ftIC  til 

Uobaooune*4 

JuatifieatlML. 


*y - - 

,  01 strltmtloa/ 
j  Availability  Co4— 

i  ]Xv»{nSS/ir 


!P\r\  j 


i 


INTRODUCTION 
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The  Institute  for  Computer  Applications  in  Science  and  Engineering  (ICASE)  is  oper¬ 
ated  at  the  Langley  Research  Center  (LaRC)  of  NASA  by  the  Universities  Space  Research 
Association  (USRA)  under  a  contract  with  the  Center.  USRA  is  a  nonprofit  consortium  of 
major  IT-  S.  colleges  and  universities. 

^  The  Institute  conducts  unclassified  basic  research  in  applied  mathematics,  numerical 
analysis,  and  computer  science  in  order  to  extend  and  improve  problem-solving  capabilities 
in  science  and  engineering,  particularly  in  aeronautics  and  space. 

ICASE  has  a  small  permanent  staff.  Research  is  conducted-pi'imarily  by  visiting  scientists 
from  universities  and  from  industry,  vrhn  hay*M‘*<d7tgnt  ^ppHr^m^nt0  for  limited  periods  of 
timej^jid-byToE^^  Members  of  NASA’s  research  staff  also  may  be  residents  at  ICASE 
m  limited  periods. 

The  major  categories  of  the  current  ICASE  research  program  are: 


Numerical  methods,  with  particular  emphasis  on  the  development  and  analysis  of  basic 
numerical  algorithms; 

Control  and  parameter  identification  problems,  with  emphasis  on  effective  numerical 
methods; 


Computational  problems  in  engineering  and  the  physical  sciences,  particularly  fluid 
dynamics,  acoustics,  and  structural  analysis; 


Computer  systems  and  software  for  parallel  computers. 


ICASE  reports  are  considered  to  be  primarily  preprints  of  manuscripts  that  have  been 
submitted  to  appropriate  research  journals  or  that  are  to  appear  in  conference  proceedings. 
A  list  of  these  reports  for  the  period  October  1,  1990  through  March  31,  1991  is  given  in  the 
Reports  and  Abstracts  section  which  follows  a  brief  description  of  the  research  in  progress. 


Presently,  ICASE  is  operated  at  NASA  Langley  Research  Center,  Hampton,  VA,  under  the  National 
Aeronautics  and  Space  Administration,  NASA  Contract  No.  NAS1-18605.  In  the  past,  support  has  been  pro¬ 
vided  by  NASA  Contract  Nos.  NAS1-18107,  NAS1-17070,  NAS1-17130,  NAS1-15810,  NAS1-16394,  NAS1- 
14101,  and  NAS1-14472. 


RESEARCH  IN  PROGRESS 


Saul  Abarbanel 

Work  continued  on  three  topics  reported  in  the  last  semiannual  report;  namely- 

•  Spurious  frequencies  due  to  numerical  boundary  treatment. 

•  Accelerating  nondiagonalizable  hyperbolic  systems  to  steady-state. 

•  Stability  of  initial-boundary  value  problems  approximated  by  high  order  Pade’-like 
finite  differences. 

In  addition  progress  has  been  made  (together  with  doctoral  student  Jeffery  Danowitz) 
on  constructing  long-wave  non-reflecting  boundary  conditions  for  the  viscous  flow  past  finite 
bodies.  At  present  the  flow  past  a  finite  flat  plate  has  been  coded  and  preliminary  results 
encourage  continuation  of  this  effort. 

Part  of  the  work  described  above  is  carried  out  together  with  David  Gottlieb  and  Mark 
Carpenter  (Fluid  Mechanics  Division,  LaRC). 


H.  T.  Banks  and  Fumio  Kojima 

Work  is  continuing  on  geometrical  heat  inverse  problems  arising  in  non-destructive  evalu¬ 
ation  methods.  Our  concern  in  this  problem  is  to  develop  an  efficient  - '  mrithm  for  identifying 
the  corrosion  profiles  :•  .*?de  a  body  using  thermographical  data.  Together  with  W.  P.  Winfree 
(Instrument  Research  Division,  LaRC),  we  are  developing  a  practical  method  for  the  estima¬ 
tion  of  the  back  surface  shape  using  the  thermal  data  from  the  front  surface.  Recently,  we 
proposed  a  spline-based  corrosion  detection  algorithm  using  the  boundary  element  method. 
Our  aim  in  this  approach  is  to  estimate  the  geometrical  structure  of  corrosions  as  well  as  to 
detect  material  flaws.  In  the  algorithm,  the  material  shape  is  characterized  by  a  periodic 
R-spline  series.  Then  the  problem  is  formulated  as  a  minimization  problem  of  the  output 
least  squares  error  subject  to  a  Volterra  integral  equation  of  the  second  kind.  Preliminary 
computational  experiments  were  successful  using  simulation  data.  We  are  currently  pursuing 
investigations  to  prove  the  existence  of  solutions  for  our  estimation  problem  and  to  analyze 
convergence  properties  for  the  proposed  computational  method. 
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H.  T.  Banks  and  Ralph  Smith 

We  are  continuing  tie  development  of  a  computationally  viable  infinite-dimensional  con¬ 
trol  and  identification  methodology  for  the  stable  and  optimal  design  of  noise  suppression 
systems.  The  basic  approach  is  to  combine  approximation  theory  with  time  domain  state 
space  modeling  to  develop  convergent  computational  algorithms  for  LQR  control  designs. 

The  physical  .  ystem  being  studied  consists  of  chambers  which  are  separated  by  a  flexible 
plate,  and  initial  efforts  have  centered  around  the  development  of  an  accurate  structure/fluid 
interaction  model  which  describes  the  transmission  of  acoustical  pressure  waves  from  one 
cylindrical  cavity  to  another  through  the  separating  flexible  plate.  This  has  led  to  the 
formulation  of  a  system  of  partial  differential  equations  consisting  of  a  3-D  wave  equation 
coupled  with  elasticity  equations  for  the  plate.  A  finite  element  scheme  is  being  developed 
to  discretize  the  system  and  the  validity  of  the  model  is  being  tested  both  numerically  and 
with  data  from  ex poiiments  by  J.  Bayer.  H.C.  Lester  and  R. J.  Silcox  (Acoustics  Division, 
LaRC). 

Once  the  modeling  and  identification  aspects  have  been  completed,  the  active  feedback 
control  of  interior  pressure  will  be  implemented  via  structural  vibration  control.  To  accom¬ 
plish  this,  feedback  to  piezo  ceramic  patches  in  the  plate  will  be  used  to  affect  plate  moments 
in  a  manner  which  optimally  controls  the  acoustic  pressure  in  the  chamber. 


Kurt  Bryan 

Work  is  underway  on  a  thermo-mechanical  inverse  problem,  namely  that  of  using  the 
mechanical  response  of  a  structure  to  external  thermal  stimuli  to  determine  its  internal 
properties  and/or  locate  defects.  The  focus  has  been  on  the  special  case  of  an  isotropic 
material  (thermally  and  mechanically)  with  an  unknown  defect  (a  hole).  Analysis  of  the 
forward  or  direct  problem  leads  to  a  system  of  partial  differential  equations  consisting  of  a 
heat  equation  coupled  to  the  linear  equations  of  elasticity.  The  thermal  source  is  assumed 
to  be  periodic  in  time,  so  that  the  elasticity  equations  also  contain  inertial  terms.  Due  to 
the  complicated  i  iture  of  the  equations,  the  computational  solution  of  the  inverse  problem 
will  likely  rely  on  an  optimization/fit-to-data  approach  requiring  an  economical  scheme  for 
rapidly  solving  the  direct  problem.  Work  so  far  has  centered  on  deriving  such  numerical 
schemes  using  the  method  of  boundary  integral  equations.  This  approach  should  be  partic¬ 
ularly  efficient  since  knowledge  of  the  mechanical  response  of  the  the  object  being  imaged 
is  required  (and  can  be  measured)  only  on  its  surface  and  not  in  its  interior.  Coding  for  a 
two-dimensional  version  of  the  problem  is  underway. 

Work  has  also  continued  on  the  inverse  conductivity  problem,  in  which  one  attempts  to 
identify  the  internal  electrical  conductivity  of  an  object  from  voltage  and  current  flux  mea- 
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surements  on  its  boundary.  The  special  case  of  identifying  a  piecewise  constant  perturbation 
of  a  constant  reference  conductivity  has  been  studied.  An  efficient  numerical  algorithm  based 
on  the  method  of  boundary  integral  equations  has  been  implemented  for  the  solution  of  the 
direct  problem,  as  well  as  a  linearized  version  of  the  direct  problem,  and  it  is  used  in  a  least 
squares  approach  for  solving  the  inverse  problem.  Performance  on  numerically  generated 
data  has  been  good  and  theoretical  convergence  results  have  been  proven. 


John  Burns 

In  earlier  work,  we  have  completed  a  study  of  several  control  problems  for  Burgers’  equa¬ 
tion.  These  nonlinear  systems  were  controlled  by  using  linear  feedback  computed  from  LQR 
theory  for  the  linearized  models.  We  established  that  under  standard  assumptions  the  finite 
element  scheme  preserved  stabilizability  and  detectability  uniformly  under  approximation. 
Numerical  results  also  indicated  that  the  LQR  controllers  smoothed  steep  gradients.  These 
results  provide  some  insight  into  the  control  of  nonlinear  distributed  parameter  systems.  We 
are  now  turning  our  attention  to  realistic  fluid  flow  problems  and  we  have  the  first  results  on 
optimal  control  of  the  full  unsteady  Navier-Stokes  flow  about  a  cylinder.  An  ICASE  report 
with  Y.  R.  Ou  is  in  preparation.  We  plan  to  continue  our  efforts  on  this  class  of  problems. 


Richard  Carter 

The  two  major  methodologies  for  improving  the  robustness  of  nonlinear  optimization 
algorithms  are  the  trust  region  approach  and  the  linesearch  approach.  Previous  work  has 
established  theory  predicting  extreme  robustness  of  trust  region  algorithms  in  the  presence 
of  noisy  function  and  gradient  evaluations.  This  robustness  has  also  been  confirmed  experi¬ 
mentally.  However,  for  the  linesearch  approach,  it  is  easy  to  generate  pathological  examples 
where  the  linesearch  algorithm  will  fail  for  even  small  amounts  of  noise  in  the  function  and 
gradient  evaluations. 

This  argues  strongly  for  the  selection  of  the  trust  region  methodology  as  the  method  of 
choice.  Although  numerical  experiments  suggest  that  the  worst-case  analysis  of  the  linesearch 
approach  is  typically  overly  pessimistic,  failures  are  still  encountered  at  levels  of  gradient 
error  for  which  the  trust  region  approach  is  still  robust. 


Leon  M.  Clancy 

The  ICASE  computing  environment  continued  to  evolve.  Two  Sun  SPARCstation  2’s 
were  purchased  along  with  6.0Gb  of  new  disk  space.  All  home  directories  formerly  on  fluke, 
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flounder,  polaris  and  croaker  were  moved  to  the  new  servers,  which  were  named  icase  and 
cobia.  The  automount  utility  was  implemented  on  all  ICASE  systems  which  support  it  in 
order  to  keep  the  home  directory  pathnames  constant  even  though  they  were  physically 
moved  to  different  fileservers. 

The  SPARCstation  2  servers  have  exceeded  expectations  as  replacements  for  the  old  Sun 
3  servers,  and  are  each  capable  of  supporting  an  additional  9.0Gb  without  stretching  the 
system’s  limits. 

The  automount  utility  implementation  mentioned  above  allowed  local  disks  scavenged 
from  ICASE’s  earliest  SPARCstation  l’s  to  be  used  as  local  swap  and  home  directories  on 
carp,  grouper,  urchin,  and  dolphin.  This  provides  a  performance  enhancement  for  those 
workstations,  yet  allows  the  users’s  whose  home  directories  reside  on  the  local  disks  trans¬ 
parent  access  to  all  ICASE  systems.  As  disks  become  available,  more  workstations  will  be 
given  this  capability. 


Thomas  W.  Crockett 

A  Sun  SPARCstation  1+  was  integrated  into  the  iPSC/860  hypercube  system  as  a  file- 
server,  compile  server,  and  remote  host.  Ihis  allows  most  program  development  activities  to 
be  offloaded  from  the  SRM,  resulting  in  improved  reliability  and  better  response  for  parallel 
jobs  running  in  the  cube.  In  addition,  a  direct  Ethernet  interface  to  the  cube  was  added,  al¬ 
lowing  network  communications  to  bypass  the  SRM,  with  performance  increased  by  a  factor 
of  5-10.  The  Ethernet  interface  also  allows  files  to  be  transferred  directly  to  and  from  the 
Concurrent  File  System. 

The  X  Window  System  and  Andrew  Toolkit  were  rebuilt  and  installed  under  SunOS  4.1.1 
on  the  ICASE  Suns. 


Thomas  W.  Crockett  and  Tobias  OrlofF 

Significant  improvements  were  made  to  the  performance  model  for  our  MIMD  render¬ 
ing  algorithm.  The  model  now  shows  good  agreement  with  experimental  results  across  a 
wide  range  of  scene  complexities  and  varying  numbers  of  processors.  Results  from  the  Intel 
iPSC/860  demonstrate  that,  while  much  of  the  communication  overhead  can  be  effectively 
overlapped  with  computation,  message  latency  is  still  the  major  overhead  on  large  num¬ 
bers  of  processors.  Hence,  distributed  memory  architectures  suitable  for  parallel  rendering 
should  (a)  allow  computation  and  communication  to  overlap,  and  (b)  have  minimal  message 
latencies. 
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Andrew  Dando 

The  behavior  of  Gortler  vortices,  and  in  particular  their  growth  rates,  have  been  inves¬ 
tigated  for  a  compressible  inviscid  flow  in  a  basic  two-dimensional  boundary  layer  in  the 
hypersonic  limit.  This  work  is  to  appear  as  a  joint  ICASE  report  with  S.  Seddougui.  The 
work  is  being  extended  to  the  case  of  a  basic  three-dimensional  boundary  layer  flow  in  an 
effort  to  determine  what  effect  crossflow  has  on  the  stability  of  these  Gortler  vortices. 


Peter  W.  Duck 

The  linear  stability  of  compressible  plane  Couette  flow  is  being  investigated.  The  proper 
basic  velocity  and  temperature  distributions  are  perturbed  by  a  small  amplitude,  normal 
mode  disturbance.  The  full  small  amplitude  disturbance  equations  are  being  solved  nu¬ 
merically  at  finite  Reynolds  numbers,  and  the  inviscid  limit  of  these  equations  is  being 
investigated  in  some  detail.  It  has  been  found  that  instability  can  occur,  although  the  sta¬ 
bility  characteristics  of  the  flow  are  quite  different  from  unbounded  flows.  The  effects  of 
viscosity  are  also  calculated,  asymptotically,  and  have  a  stabilizing  role  in  all  the  cases  inves¬ 
tigated.  Exceptional  regimes  to  the  problem  occur  when  the  wavespeed  of  the  disturbances 
approaches  the  velocity  of  either  of  the  walls,  and  these  regimes  are  also  being  analyzed  in 
some  detail.  Finally,  the  effect  of  imposing  radiation-type  boundary  conditions  on  the  up¬ 
per  (moving)  wall  (in  place  of  impermeability)  is  being  investigated,  and  this  yields  results 
common  to  both  bounded  and  unbounded  flows.  The  research  is  being  done  in  collaboration 
with  G.  Erlebacher  and  M.  Y.  Hussaini. 


Thomas  Eidson 

A  project  to  study  the  performance  of  typical  CFD  codes  on  the  Intel  iPSC/860  was 
begun.  The  study  includes  single  node  performance  as  well  as  the  modification  of  typical 
CFD  algorithms  for  parallel  performance  (better  load  balance  and  communication).  In 
collaboration  with  Joel  Saltz  and  Kay  Crowley,  the  code  transformations  determined  useful 
in  the  above  study  are  being  used  as  design  goals  to  evaluate  and  possibly  develop  code 
transformation  tools  and  compilers. 

Initially,  the  work  focused  on  studying  basic  hypercube  operations-  node  communication 
and  single  node  performance  (cache  management  and  assembly  coding  of  vector  operations). 
A  vector  library  was  developed  (partially)  and  a  programming  strategy  for  that  library  was 
studied.  When  the  relevant  data  is  in  the  data  cache,  the  vector  library  routines  can  attain  37 
Mfiops  (40  theoretical  peak)  for  simple  vector  operations  and  55  Mflops  (80  peak)  for  chaining 
operations.  In  more  realistic  kernels  with  typical  data  flow  patterns  where  sufficient  data 
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does  not  always  stay  in  cache,  rates  vary  from  7  to  20  Mflops.  Node  communication  will 
reduce  these  single  node  rates  further  for  a  typical  code.  These  single  node  rates  using  the 
library  are  faster  than  those  attained  with  the  PGI  compiler  but  the  speed-up  is  usually  less 
than  2.  The  PGI  compiler  does  partial  vectorization  and  has  a  limited  vector  library.  Work 
is  continuing  on  optimizing  the  current  vector  library/programming  strategy  and  studying 
other  approaches. 

An  explicit,  finite-difference  code  developed  by  Gordon  Erlebacher  to  study  3-D,  com¬ 
pressible  turbulence  is  being  converted  to  the  hypercube  to  evaluated  various  performance 
programming  ideas.  The  first  kernel  being  studied  is  a  periodic  tridiagonal  solver  with  mul¬ 
tiple  right-hand  sides.  Several  data  structures  and  algorithm  approaches  are  being  coded 
and  studied. 


Gordon  Erlebacher 

Work  on  shear  flow  compressible  turbulence  in  collaboration  with  S.  Sarkar  and  M.Y. 
Hussaini  has  been  completed.  The  shear  flow  code  has  been  debugged  and  first  results 
analyzed.  A  report  has  been  written  and  accepted  in  the  journal  of  Theoretical  and  Compu¬ 
tational  Fluid  Dynamics.  Eigenvalue/eigenvector  analysis  combined  with  statistical  methods 
is  used  to  analyze  the  shear  flow  data.  The  objective  is  to  correlate  the  visual  pictures  of 
streaks  with  more  quantitative  diagnostics. 

In  collaboration  with  M.Y.  Hussaini  and  M.  Malik  and  C.-L.  Chang  (High  Technology 
Inc.),  the  Parabolized  Stability  formulation  applied  to  supersonic  flows  over  boundary-layers 
has  been  completed.  Both  the  linear  and  the  nonlinear  PSE  has  been  implemented  and 
checked.  The  nonlinear  2-D  PSE  was  checked  against  second  mode  data  generated  in  1988. 
Results  of  this  work  will  be  presented  in  Hawaii.  The  programming  wss  performed  by  C.-L. 
Chang. 

Together  with  T.A.  Zang  (Fluid  Dynamics  Division,  LaRC)  and  D.  Pruett  (National 
Research  Council  Fellow),  direct  numerical  simulations  of  supersonic  transitional  flow  over 
boundary-layers  were  completed.  Simulations  were  performed  of  flows  over  sharp  cones  at 
zero  angle  of  attack  at  a  free-stream  Mach  number  of  8.  Shape  factors  were  computed  and 
qualitative  comparisons  of  Reynolds  stress  plots  with  experimental  results  were  done.  Good 
qualitative  agreement  was  obtained.  An  ICASE  report  was  completed. 

Work  with  P.  Duck  and  M.Y.  Hussaini  on  the  stability  of  compressible  plane  Couette  flow 
has  produced  new  types  of  instabilities.  This  work  combines  asymptotics  of  high  Reynolds 
number  flows  with  linear  stability  calculations  of  the  linearized  Navier-Stokes  equations  as 
finite  Reynolds  number  to  clarify  some  of  the  mechanisms  that  lead  to  unstable  modes  in 
the  supersonic  regimes.  In  this  work,  we  tried  to  relate  the  modes  found  in  Couette  flow 
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with  those  known  to  exist  in  supc  sonic  boundary  layers.  An  ICASB  report  describing  this 
work  has  been  completed. 

Work  is  in  progress  with  S.  Biringen  at  Colorado  Stele  University  on  the  Direct  Numerical 
Simulation  of  compressible  Couette  flow. 


James  F.  Geer 

Work  is  continuing  on  hybrid  perturbation/ variational  techniques  for  solving  a  variety 
of  problems  for  both  ordinary  and  partial  differential  equations.  For  ODE’s,  we  are  cur¬ 
rently  applying  the  methods  to  several  classes  of  nonlinear  vibration  problems,  with  special 
emphasis  on  resonant  frequency  calculations  for  systems  of  equations.  For  PDE’s,  we  are 
investigating  how  the  method  might  be  applied  to  some  exterior  boundary  value  problems 
for  elliptic  PDE’s,  when  the  boundary  of  the  domain  D  has  an  irregular  shape.  The  basic 
idea  is  to  treat  D  as  a  perturbation  of  a  simple  domain  R  (e.g.  a  circle  or  an  ellipse,  in  two 
dimensions)  by  embedding  it  in  a  one  parameter  family  of  domains  indexed  by  a  parameter 
a,  where  0  <  a  <  1.  Currently,  we  are  applying  the  method  to  some  eigenvalue  problems  for 
irregularly  shaped  domains  and  to  some  problems  of  flows  about  geometrically  complicated 
bodies.  In  addition,  we  are  combining  our  method  with  some  homotopy  methods  to  further 
enhance  the  accuracy  of  the  results  we  obtain.  In  particular,  we  hope  this  will  lead  to  a 
new  extension  of  slender  body  theory.  So  far,  the  preliminary  results  have  been  done  with 
C.  Andersen  (College  of  William  and  Mary),  while  possible  applications  are  being  discussed 
with  E.  Liu  (Fluid  Mechanics  Division,  LaRC)  and  M.  Hemsch  (Lockheed  Engineering  and 
Science  Co.). 

Investigations  are  also  continuing  concerning  some  fundamental  properties  of  a  class  of 
“almost”  singular  integral  equations  of  the  first  find  which  are  useful  in  representing  solutions 
to  certain  elliptic  exterior  boundary  value  problems.  These  equations  typically  have  the 
property  that  the  domain  of  integration  R  is  a  proper  subset  of  the  domain  of  validity  D 
of  the  equation.  Special  consideration  is  being  given  to  the  idea  of  analytically  continuing 
the  solution  into  the  domain  D.  In  fact,  this  has  now  been  done  for  a  large  class  of  one¬ 
dimensional  integrals,  such  els  those  which  occur  in  the  representation  of  solutions  involving  a 
body  of  revolution.  Application  of  the  results  to  several  two  and  three  dimensional  problems 
involving  slender  or  thin  bodies  are  being  carried  out.  The  symbolic  manipulation  systems 
MACSYMA  and  Mathematica  have  been  used  in  some  of  the  preliminary  investigations. 

The  problem  of  describing  aerodynamically  generated  sound  from  compact  sources  of 
vorticity  is  being  studied  from  a  perturbation  point  of  view,  with  the  eddy  Mach  number 
appearing  els  the  perturbation  parameter.  The  problem  appears  to  be  well  suited  to  an  ap¬ 
plication  of  a  slightly  modified  version  of  the  method  of  multiple  scales.  In  particular,  the 
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problem  of  determining  the  sound  generated  by  a  sphere  executing  arbitrary  (small)  oscilla¬ 
tions  in  zero  mean  flow  conditions  has  been  investigated  successfully  by  this  technique.  Once 
the  perturbation  solution  is  fully  understood,  the  hybrid  perturbation  Galerkin  technique 
will  be  applied  to  it,  with  the  goal  of  extending  the  usefulness  (and  accuracy)  of  the  solution 
to  higher  Mach  numbers  than  is  possible  using  the  perturbation  solution  alone.  This  work 
is  being  carried  out  with  J.  Hardin  (Acoustics  Division,  LaRC). 


David  Gottlieb 

The  issue  of  asymptotic  stability,  as  opposed  to  Lax  stability,  becomes  more  and  more 
important  for  long  time  simulations.  The  notion  of  Lax  stability  permits  a  stable  scheme  to 
have  a  solution  that  grows  exponentially  in  time;  thus  a  Lax  stable  scheme  may  be  useless 
for  the  numerical  simulations  of  genuinely  time  dependent  problems.  We  have  found  that  for 
high  order  schemes,  in  particular  fourth  and  sixth  order  compact  finite  difference  schemes, 
the  numerical  inflow  boundary  condition  often  induces  time  growth  which  renders  long  time 
simulations  useless.  We  have  developed  a  theory  for  checking  asymptotic  stability  for  systems 
of  hyperbolic  equations  and  used  it  to  verify  that  a  new  boundary  treatment  based  on  the 
penalty  method  is  stable.  This  work  is  done  with  D.  Carpenter  (Fluid  Mechanics  Division, 
LaRC)  and  S.  Abarbanel. 

We  have  continued  our  work  on  spectral  methods  for  shock  wave  calculations.In  this 
context  we  found  a  method  that  completely  removes  the  Gibbs  phenomenon,  even  at  the 
discontinuity.  Numerical  simulations  of  two  dimensional  interactions  of  shocks  and  vortices 
demonstrate  the  applicability  of  spectral  methods  for  those  problems. 

We  are  looking  at  domain  decomposition  methods  as  means  to  parallelize  spectral  meth¬ 
ods.  To  minimize  communications  we  use  explicit  methods  to  advance  the  interface  points 
and  implicit  methods  for  the  inner  domains.  A  substantial  saving  has  been  observed  for 
some  parabolic  problems. 


Philip  Hall 

Hypersonic  boundary  layers  can  be  unstable  for  a  variety  of  different  instability  mecha¬ 
nisms.  We  have  been  investigating  the  interaction  of  inviscid,  Gortler  and  Tollmien  Schlicht- 
ing  waves  in  Sutherland  Law  fluids.  The  effect  of  real  gas  properties  and  shocks  on  the 
instabilities  has  been  studied,  as  has  the  effect  of  compliant  walls  on  vortex  instabilities.  A 
new  class  of  exact  Navier  Stokes  solutions  has  been  derived  and  we  have  shown  the  relevance 
of  the  solutions  to  large  amplitude  instabilities  in  boundary  layer  flows.  Other  work  in  three 
dimensional  boundary  layer  instability  theory  concerns  vortex  wave  interactions  and  the 
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interaction  of  crossflow,  Gortler  and  Tollmien  Schlichting  modes.  The  receptivity  problem 
for  nonparallel  modes  in  disc  flows  has  been  solved.  Work  is  also  being  carried  out  on  the 
secondary  instability  of  stronglj  nonlinear  vortex  flows  in  incompressible  and  compressible 
boundary  layers.  This  work  has  already  quantitatively  reproduced  many  of  the  key  features 
of  the  breakdown  structure  of  streamwise  vortices;  present  work  is  concentrated  on  the  hy¬ 
personic  limit  of  that  problem.  Work  on  the  instability  of  heated  boundary  layers  is  also  in 
progress. 


Thomas  L.  Jackson 

Work  focuses  on  the  time  evolution  on  the  interaction  of  a  diffusion  flame  with  a  vortex, 
which  directly  applicable  to  the  large  scale  structures  found  in  compressible  mixing  layers.  In 
particular,  the  ignition  process  and  the  resulting  structure  are  under  current  investigation. 
A  combination  of  asymptotics  and  numerics  is  used  to  reduce  the  complex  problem  to  a 
model  problem,  thus  isolating  key  physical  effects  for  analysis.  This  work  is  in  collaboration 
with  M.  Macaraeg  (Fluid  Mechanics  Division,  LaRC)  and  M.Y.  Hussaini. 


Peter  A.  Jacobs 

Work  has  continued  on  the  simulation  of  transient,  hypervelocity  flows  in  shock  tunnels. 
We  attempted  to  simulate  the  starting  flow  in  a  Mach  8  axisymmetric  nozzle  and  reported 
the  results  in  ICASE  Report  91-1.  This  attempt  was  only  partially  successful  and  we  are 
currently  trying  again  with  an  finite- volume  upwinding  code  (see  Jacobs  and  Scroggs).  De¬ 
spite  a  number  of  deficiencies,  the  simulations  proved  to  be  useful  in  guiding  the  analysis  of 
experimental  measurements.  Two  sets  of  previously  obtained  measurements  have  since  been 
reexamined.  ICASE  Report  91-24  documents  some  of  the  experimental  data  for  a  Mach  4 
nozzle  and  a  larger  Mach  8  nozzle  while  ICASE  Interim  Report  16  documents  the  data  for 
a  scramjet  combustor  experiment. 


Peter  A.  Jacobs  and  Jeffrey  S.  Scroggs 

We  have  developed  yet  another  program  for  the  integration  of  the  Navier-Stokes  equations 
on  a  two-dimensional  structured  mesh.  Although  the  code  currently  considers  only  a  single¬ 
block  grid,  the  full  version  of  the  code  will  describe  a  flow  domain  as  a  set  of  abutting  blocks, 
each  consisting  of  a  tensor-product  mesh  of  quadrilateral  cells.  The  flow  field  is  recorded  as 
cell-average  values  at  cell  centers  and  explicit  time-stepping  is  used  to  update  conserved 
quantities.  Second-order  flux-difference  splitting  (via  an  Osher-type  Riemann  solver  and 
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limiter)  is  used  to  calculate  inviscid  fluxes  across  cell  faces  while  central  differences  (via  the 
divergence  theorem)  are  used  to  calculate  the  viscous  fluxes.  The  unusual  features  of  this 
program  are  that  it  is  written  in  C  and  makes  extensive  use  of  sophisticated  data  structures 
to  encapsulate  the  data.  The  idea  of  writing  the  code  this  way  is  to  make  it  easier  (than 
traditional  FORTRAN  codes)  to  “parallelize”  for  the  MIMD  style  of  parallel  computer.  The 
code  is  intended  to  be  a  test-bed  for  implementation  on  parallel  computers  and  work,  in 
collaboration  with  J.  Saltz,  will  continue  in  the  later  half  of  the  year. 


Harry  F.  Jordan 

Multiprocessors  are  currently  manifest  in  two  distinct  forms:  shared  memory  and  dis¬ 
tributed  memory.  The  distinguishing  feature  is  the  level  at  which  system  latency  resulting 
from  data  movement  is  handled.  We  are  studying  the  software  support  for  suppressing 
the  distinctions  among  data  based  on  its  location  in  a  message  passing  multiprocessor  and 
relating  the  software  techniques  for  algorithm  level  global  referencing  to  the  hardware  tech¬ 
niques  used  to  support  a  shared  address  space  at  the  machine  instruction  set  level.  The  key 
issue  is  the  management  of  latency  involved  in  data  movement.  Latency  can  be  reduced 
by  appropriate  placement  of  data  or  masked  by  overlapping  useful  computation  with  data 
movement.  Data  reduction  is  done  by  problem  mapping  in  software  for  distributed  memory 
multiprocessors  and  by  caches  in  hardware  for  shared  memory  machines.  Latency  masking  is 
done  by  scheduling  message  traffic  in  advance  of  the  need  for  its  data  in  distributed  memory 
machines  and  by  various  forms  of  pipelining  in  the  hardware  of  shared  memory  multipro¬ 
cessors.  Comparing  and  contrasting  the  techniques  used  for  latency  management  in  both 
domains  can  yield  new  insight  into  the  best  combination  of  both  levels  of  support  for  large 
scale  multiprocessors. 

In  a  separate  effort,  we  are  studying  speed  scalable  optical  computer  architectures.  These 
architectures  postulate  that  the  speed  of  signal  propagation  between  gates  is  commensurate 
with  logic  switching  speeds.  Since  the  speed  of  light  is  the  upper  limit  on  information 
transfer  speed,  optical  computing  forms  a  natural  domain  to  study  the  fine  grain  pipelined 
architectures  which  result  from  this  postulate.  An  essential  element  o.  these  speed  of  light 
architectures  is  relaxing  the  assumption  of  instantaneous  signal  propagation  within  a  sequen¬ 
tial  circuit  clock  cycle.  We  are  both  developing  algorithms  based  on  the  above  assumptions 
and  working  with  prototype  hardware  to  explore  this  new  design  domain.  The  connection 
between  this  work  and  that  mentioned  above  is  the  acceptance  of  latency  as  an  inherent 
feature  of  very  large  or  very  fast  systems  and  the  use  of  creative  methods  to  manage  that 
latency. 
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David  E.  Keyes 

The  divide-and-conquer  paradigm  of  iterative  domain  decomposition,  or  substructuring, 
has  become  a  practical  tool  in  computational  fluid  dynamics  applications  because  of  its  flex¬ 
ibility  in  accommodating  adaptive  refinement  through  locally  uniform  (or  quasi-uniform) 
grids,  its  ability  to  exploit  multiple  discretizations  of  the  operator  equations,  and  the  mod¬ 
ular  pathway  it  provides  towards  parallelism.  Our  recent  work  at  ICASE  (summarized  in 
ICASE  Reports  91-19  and  91-20  and  related  references  therein)  has  focused  on  a  class  of 
domain-decomposed  preconditioners  based  on  non-overlapping  subregions  and  a  coarse  grid 
of  subdomain  vertices.  We  have  evaluated  the  algebraic  convergence  rate  and  the  paralleliz- 
ability  of  algorithms  based  on  such  preconditioners  for  nonlinear  multicomponent  systems 
arising  from  fluid  transport  applications  and  have  demonstrated  how  various  local  refine¬ 
ments  can  be  incorporated  into  the  composite  discrete  operator.  The  parallel  testbeds  of 
primary  interest  are  the  medium  scale  MIMD  machines,  such  as  the  Intel  iPSC/860. 


D.  Glenn  Lasseigne 

Our  research  focuses  on  the  interactions  of  disturbances  in  the  flow  field  with  both  react¬ 
ing  and  non-reacting  shocks.  A  combination  of  asymptotics  and  numerics  is  used  to  reduce 
complex  problems  to  model  problems,  thus  isolating  key  physical  effects  for  analysis.  Spe¬ 
cific  problems  include:  the  non-linear  interactions  of  vorticity /detonation  waves,  the  effects 
of  heat  release  on  the  response  of  an  oblique  detonation  subjected  to  wedge  oscillations,  and 
the  coupling  of  upstream  disturbances  and  wedge  oscillations  on  the  stability  of  an  oblique 
detonation.  The  goal  is  to  develop  numerical  routines  which  capture  the  instability  modes 
predicted  by  linear  theory.  Other  problems  include  the  effects  of  streamwise  vorticity  on  the 
stagnation  point  flow  and  the  effects  of  more  general  disturbances  on  a  flame  in  the  viscous 
stagnation  region. 


Dimitri  Mavriplis 

Work  is  continuing  on  the  use  of  unstructured  meshes  for  solving  computational  fluid 
dynamics  problems  in  both  two  and  three  dimensions. 

In  two  dimensions,  a  viscous  turbulent  flow  solver  has  been  developed.  Work  has  con¬ 
centrated  on  the  implementation  and  modification  of  a  multiple  field  equation  turbulence 
model  into  the  two-dimensional  code  (ICASE  Report  91-11).  Research  on  implicit  methods 
for  two-dimensional  viscous  turbulent  flows  has  also  been  pursued.  An  ILU  (incomplete  LU 
factorization)  preconditioned  GMRES  implicit  iterative  solver  has  been  developed  in  con¬ 
junction  with  V.  Venkatakrishnan  (Computer  Science  Corporation,  Moffett  Field,  CA)  for 
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the  two-dimensional  Euler  and  turbulent  Navier-Stokes  equations  on  unstructured  meshes. 

In  three  dimensions,  an  inviscid  flow  solver  has  been  developed  incorporating  adaptive 
meshing  techniques  and  an  unstructured  multigrid  algorithm,  thus  enabling  the  efficient  and 
accurate  computation  of  three-dimensional  flow  fields  about  complex  configurations.  Future 
work  will  be  directed  at  implementing  the  viscous  turbulent  prediction  capability  developed 
in  the  two-dimensional  framework  into  three  dimensions. 

In  conjunction  with  J.  Saltz,  S.  Das,  and  R.  Ponnusamy,  the  single-grid  non-adaptive 
three-dimensional  unstructured  Euler  solver  has  been  implemented  on  the  Intel  iPSC/860 
using  the  PARTI  primitives.  This  work  is  aimed  at  demonstrating  the  applicability  of  the 
PARTI  primitives  to  practical  implementations,  and  also  to  employ  experience  gained  from 
particular  implementations  to  aid  in  the  design  or  modification  of  the  primitives.  The  im¬ 
plementation  will  be  extended  in  the  future  to  include  the  multigrid  acceleration  capability 
and  adaptive  meshing  strategies. 


Piyush  Mehrotra 

Current  programming  environments  for  distributed  memory  machines  force  the  program¬ 
mer  to  decompose  large  data  structures  into  a  collection  of  pieces,  each  piece  being  “owned” 
by  a  single  process.  The  interaction  between  pieces  is  then  explicitly  specified  using  the 
underlying  message  passing  primitives.  Decomposing  data  structures  in  this  way  and  explic¬ 
itly  specifying  communication  not  only  is  complex  and  error  prone  but  also  “hardwires”  the 
distributions  and  algorithmic  choices. 

We  have  been  developing  a  high-level  programming  environment,  called  Kali,  to  alleviate 
some  of  these  problems.  In  our  approach,  the  code  is  specified  using  a  global  name  space  in 
a  distribution  independent  manner.  The  user,  however,  controls  performance  critical  issues 
such  as  load  balancing  by  specifying  the  distribution  of  data  through  high-level  annotations. 
We  have  run  experiments  to  show  that  by  using  our  approach,  it  is  easy  for  the  user  to 
experiment  with  alternative  data  distributions  with  minimal  changes  to  the  source  code. 
This  is  because  the  complex  details  of  interprocessor  communication  is  the  responsibility 
of  the  compiler  and  the  run-time  environment  rather  than  the  programmer.  Our  system 
has  provision  for  both  static  and  dynamic  distributions  and  we  are  continuing  to  investigate 
transformations  and  optimizations  needed  to  efficiently  support  both  within  the  framework 
of  our  system. 
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Kirsten  Morris 

Uniform  stabilizability  of  approximations  appears  to  be  critical  to  justifying  the  use 
of  finite-dimensional  approximations  in  controller  design  for  infinite-dimensional  systems. 
Several  general  results  exist  for  systems  with  bounded  input  operators.  However,  in  boundary 
control  of  partial  differential  equation  systems,  the  input  operator  is  unbounded.  A  study  of 
the  requirements  on  the  approximation  scheme  in  order  for  uniform  stabilizability  to  hold  is 
being  done. 

Accelerometers  are  gaining  increasing  popularity  as  sensors  in  control  of  structures.  They 
are  low  in  cost,  do  not  need  to  be  frequently  recalibrated  and  are  accurate.  They  are  of 
particular  importance  in  control  of  large  space  structures,  since  unlike  displacement  sensors, 
no  frame  of  reference  is  needed.  Despite  their  practical  importance,  the  theoretical  aspects 
*  of  accelerometers  have  not  been  explored. 

“Well-posed”  control  systems  have  several  desirable  characteristics:  (1)  the  map  from 
input  to  state  is  well-defined  and  continuous  in  some  sense  and,  (2)  the  map  from  the  state 
to  the  output  is  also  well-defined  and  continuous.  In  addition,  an  equivalence  exists  between 
internal  exponential  stability  and  external  stability  for  stabilizable/detectable  well-posed 
systems.  It  can  be  shown  that  structural  control  systems  with  accelerometers,  are  not  in 
general  well-posed  in  the  sense  of  Salamon.  This  is  not  surprising,  since  the  output  is  not  a 
direct  measurement  of  a  state. 

However,  an  input-output  relation  and  an  equivalence  between  time  and  frequency  do¬ 
main  behavior  of  these  systems  can  be  obtained.  This  result  justifies  the  use  of  frequency 
domain  methods  in  the  design  of  controllers  for  such  systems.  Furthermore,  approximation 
methods  based  on  distance  in  a  graph  metric  may  be  used  in  controller  design.  This  analysis 
also  provides  insight  into  the  nature  of  the  “generalized  transfer  function”  which  exists  for 
well-posed  systems.  Work  with  H.T.  Banks  to  determine  whether  an  equivalence  between 
internal  and  external  stability  exists  for  these  control  systems  is  underway. 

Collaboration  with  J.  Juang  (Structures  and  Dynamics  Division,  LaRC)  is  directed  at 
design  techniques  for  structural  control  which  do  not  rely  on  an  accurate  damping  model. 
Due  to  the  practical  advantages  of  accelerometers,  we  are  currently  focusing  our  efforts  on 
structures  with  acceleration  feedback.  J.  Juang  and  A.  Bruner  (Lockheed  Engineering  and 
Science  Co.)  have  implemented  an  unconditionally  stable  controller  which  improves  the 
closed  loop  damping.  We  are  currently  attempting  to  generalize  this  approach  to  obtain 
a  closed  loop  system  which  possesses  a  small  settling  time  over  a  wide  range  of  open  loop 
structural  damping  behavior. 
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Naomi  Naik  and  John  Van  Rosendale 

Multigrid  algorithms  can  solve  an  elliptic  PDE  on  an  n- point  mesh  in  0(n)  sequential 
operations,  making  them  optimal  order  sequential  algorithms.  They  are  also  effective  parallel 
algorithms;  V-cycle  multigrid  can  solve  elliptic  PDEs  on  an  n-point  mesh  in  0(log2  n)  parallel 
steps.  However,  though  this  parallel  complexity  is  reasonable,  it  is  a  factor  of  log  n  worse  than 
optimal.  This  is  a  direct  consequence  of  the  “idle  processor  problem”:  parallel  architectures 
do  not  have  enough  useful  work  to  do  on  coarse  grid  levels. 

The  Frederickson-McBryan  Parallel  Multigrid  algorithm  addresses  the  “idle  processor 
problem,”  and  improves  the  constant  in  the  0(log2  n)  complexity  bound  somewhat,  but 
does  not  change  the  order  of  complexity.  The  only  apparent  way  of  changing  this  is  to  use 
concurrent  iteration,  in  which  all  multigrid  levels  are  relaxed  simultaneously. 

In  recent  work  we  have  combined  the  idea  of  concurrent  iteration  with  use  of  Fredrickson- 
McBryan-style  multiple  coarse  grids,  to  get  a  concurrent  relaxation  algorithm  which  is 
relatively  easy  to  analyze.  The  resulting  algorithm  has  parallel  complexity  O(log“n),  for 
a  G  (1,  2),  where  the  specific  value  of  a  depends  on  the  interpolation  and  smoothing  oper¬ 
ators  used.  This  is  the  best  parallel  complexity  bound  for  multigrid  obtained  so  far. 

In  another  recent  development,  we  have  constructed  a  “robust”  point  relaxation  multigrid 
algorithm  using  different  kinds  of  multiple  coarse  grids.  As  is  well  known,  multigrid  algo¬ 
rithms  based  on  line  or  plane  relaxation  in  each  of  the  coordinate  directions  can  be  “robust,” 
in  the  sense  that  their  convergence  rate  is  independent  of  mesh  aspect  ratio.  However,  such 
algorithms  are  expensive,  and  parallelize  poorly.  Consider,  for  example  a  three  dimensional 
Chebyshev  grid.  In  this  case,  only  plane-relaxation  multigrid  suffices.  Such  algorithms  have 
limited  and  awkward  parallelism,  requiring  0( log3  n)  time  per  V-cycle,  or  at  least  0(log4  n) 
time  to  convergence. 

An  alternative  way  of  achieving  robustness,  which  avoids  line  and  plane  relaxations, 
is  to  use  different  kinds  of  multiple  coarse  grids,  formed  by  semi-coarsening  separately  in 
each  of  the  coordinate  directions.  This  is  an  idea  due  to  Wim  Mulder,  who  studied  its 
use  in  the  context  of  a  grid  alignment  problem  for  hyperbolic  PDEs.  We  discovered  a 
simple  modification  of  Mulder’s  algorithm  for  the  elliptic  case,  which  leads  to  a  fast  and 
effective  robust  multigrid  algorithm.  The  new  algorithm  is  competitive  with  line  relaxation 
on  sequential  architectures,  and  far  more  attractive  on  parallel  machines. 

One  interesting  idea  here  is  to  combine  this  algorithm  with  concurrent  iteration.  The 
resulting  algorithms  map  very  well  to  a  variety  of  parallel  architectures,  including  SIMD 
machines  like  the  Connection  Machine.  In  joint  work  with  Joe  Dendy  of  Los  Alamos  National 
Laboratory,  we  are  studying  the  effectiveness  of  this  approach  on  a  variety  of  difficult  elliptic 
problems,  including  anisotropic  problems  with  large  coefficient  jumps.  We  are  also  exploring 
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the  generalization  of  this  and  other  concurrent  iteration  multigrid  algorithms  to  nonlinear 
PDEs.  A  concurrent  iteration  FAS-style  algorithm  has  been  defined,  and  will  soon  be  tried 
on  a  sequence  of  nonlinear  test  problems. 


David  Nicol 

Together  with  A.  Rifkin  (College  of  William  and  Mary)  we  have  developed  and  are 
studying  an  algorithm  for  optimally  partitioning  a  two  dimensional  workload  in  a  rectangle 
into  rectangles  such  that  each  rectangle  has  only  one  neighboring  rectangle  in  each  of  the 
four  directions  (NEWS)  of  the  compass.  This  property  is  important  on  parallel  architectures 
where  the  cost  of  communicating  with  a  non-NEWS  neighbor  is  extremely  high. 

Algorithms  for  executing  trace-driven  cache  simulations  on  SIMD  architectures  are  un¬ 
der  study  with  A.  Greenberg  and  B.  Lubachevsky  (AT&  T,  Bell  Laboratories).  We  have 
developed  algorithms  suitable  for  simulating  a  general  class  of  stack  replacement  policy  al¬ 
gorithms,  as  well  as  random  replacement. 

Work  on  load  balancing  is  continuing  with  R.  Simha  (College  of  William  and  Mary)  and 
D.  Towsky  (University  of  Massachusetts).  We  are  looking  at  the  application  of  the  theory 
of  majorization  to  problems  in  statically  load  balancing  a  stochastically  evolving  workload. 
Majorization  permits  us  to  assert  that  one  mapping  is  better  than  another  in  a  variety  of 
ways,  including  expected  execution  time,  variation  in  execution  time,  space-time  product, 
and  reliability. 

In  a  separate  effort  with  D.  Reed,  we  are  developing  simulation  models  and  analytic 
models  of  dynamic  load-balancing  policies  that  choose  a  destination  processor  at  random. 
Such  policies  are  easy  to  implement,  although  they  ignore  state  information.  Our  aim  is  to 
assess  the  quality  and  costs  of  these  types  of  policies. 


Yuh-Roung  Ou 

In  collaboration  with  John  Burns  work  is  continuing  on  the  development  of  theoretical 
methods  and  numerical  algorithms  for  various  optimal  boundary  control  problems  arising  in 
fluid  flow.  The  main  thrust  of  our  current  investigation  is  on  simulation  and  control  of  an 
unsteady  flow  generated  by  a  circular  cylinder  undergoing  a  combined  (steady  or  unsteady) 
rotatory  and  rectilinear  motion. 

By  treating  the  rotation  rate  as  a  control  variable  in  this  model,  the  fundamental  problems 
associated  with  the  development  of  the  alternate  shedding  of  vortices  and  the  induced  body 
forces  on  the  cylinder  axe  studied  numerically.  Although  this  research  will  cover  a  wide  range 
of  time-dependent  rotation  rates,  the  computational  experiments  so  far  have  been  performed 
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under  the  following  three  kinds  of  rotation:  (1)  constant  speed  of  rotation,  (2)  time-harmonic 
rotatory  oscillation,  and  (3)  time-periodic  rotation. 

We  have  tested  the  model  against  the  experimental  work  of  Coutanceau  &  Menard  with 
excellent  agreement.  The  numerical  results  elucidate  the  significant  influence  of  time- varying 
rotation  upon  the  characteristics  of  wake  and  the  process  of  vortex  shedding.  Precise  un¬ 
derstanding  of  this  mechanism  in  boundary  layer  control  may  provide  an  effective  way  for 
lift  enhancement  and  reduction  of  vortex-induced  oscillation.  The  goal  of  this  work  is  to 
contribute  a  useful  systematic  algorithm  for  many  practical  control  design  problems  of  high- 
performance  aero/hydro-vehicles,  and  to  lay  a  bridge  between  control  theory  and  the  control 
of  fluid  flow  systems. 


Demetrious  T.  Papageorgiou 

Work  has  been  continuing  on  the  stability  of  compressible  wakes  and  shear  layers.  In 
particular  the  hypersonic  limit  of  the  wake  stability  problem  was  analyzed,  the  main  finding 
being  that  the  dominant  characteristics  are  derivable  from  a  local  analysis  of  the  vorticity 
mode  in  the  vicinity  of  the  generalized  inflection  point  of  the  basic  flow.  A  novel  hypersonic 
stability  equation  was  also  derived  and  analyzed,  valid  in  the  far-wake  limit  where  physically, 
the  vorticity  layer  is  no  longer  asymptotically  thin.  Other  on-going  research  includes  the 
stability  of  shear  layers  with  a  wake  component  (with  M.Macaraeg,  Fluid  Mechanics  Division, 
LaRC)  and  asymptotic  long- wave  stability  of  compressible  shear  layers  (with  T.  Jackson). 

Work  has  also  been  completed  (with  Y.  Smyrlis)  to  characterize  transitions  to  chaos  for 
the  Kuramoto-Sivashinsky  equation,  a  dissipative  partial  differential  equation.  We  managed 
to  accurately  compute  classical  routes  to  chaos  according  to  the  theory  of  period-doubling  bi¬ 
furcations  and  have  found  excellent  agreement  with  the  dynamics  of  one  dimensional  discrete 
maps;  in  particular,  we  have  computed  both  of  Feigenbaum’s  universal  constants.  Work  is 
under  way  on  two  fronts  (i)  the  inclusion  of  dispersive  effects  in  the  equation,  (ii)  the  predic¬ 
tion  of  the  dynamics  at  different  physical  parameter  values  from  those  in  the  period-doubling 
cascade. 


Ugo  Piomelli 

Initial  studies  for  the  application  of  a  new  model  in  large-eddy  simulation  of  compressible 
flows  have  been  performed.  The  dynamic  eddy  viscosity  model  of  Germano,  Piomelli,  Moin 
and  Cabot  [ Phys .  Fluids  A ,  3(7),  1991]  has  been  reformulated  for  compressible  flows,  and  its 
extension  to  the  energy  equation  has  been  derived.  This  model  utilizes  the  energy  content  of 
the  smallest  resolved  scales  to  determine  the  model  constants  (analogous  to  the  Smagorinsky 
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constant  and  to  the  turbulent  Prandtl  number)  dynamically.  The  model  constants  are,  thus, 
capable  of  adjusting  to  the  local  state  of  the  flow.  This  feature  has  been  successfully  exploited 
for  the  simulation  of  both  turbulent  and  transitional  flows.  A  priori  tests  and  simulations 
of  compressible  homogeneous  isotropic  turbulence  are  under  way. 


Peter  Protzel 

In  a  joint  effort  with  D.  Palumbo  (Information  Systems  Division,  LaRC),  we  successfully 
completed  an  initial  study  of  the  fault-tolerance  of  a  certain  type  of  Artificial  Neural  Net¬ 
works  (ANNs)  that  can  be  used  to  solve  optimization  problems.  One  application  that  we 
investigated  is  an  ANN  that  controls  the  reallocation  of  tasks  in  a  fault-tolerant,  distributed 
multiprocessor  system.  This  reallocation  is  necessary  if  one  of  the  processors  fails  and  all 
tasks  have  to  be  distributed  among  the  remaining  processors  as  quickly  as  possible  by  ob¬ 
serving  certain  constraints  and  by  approximately  balancing  the  load  of  the  processors.  Thus, 
the  task  allocation  is  a  constrained  optimization  problem  that  can  be  solved  by  an  ANN. 
However,  the  ANN  then  becomes  a  critical  component  of  the  system  and  the  fault-tolerance 
of  the  ANN  itself  is  of  major  concern.  In  order  to  study  the  fault-  tolerance,  we  simulated 
the  operation  of  the  ANN  in  the  presence  of  different  types  of  component  failures.  Our 
results  show  a  surprising  degree  of  fault-tolerance  with  only  a  slight  performance  degrada¬ 
tion  even  after  multiple  faults.  A  discussion  of  the  results  can  be  found  in  an  upcoming 
NASA  Technical  Paper  and  an  ICASE  Report  currently  in  preparation.  The  conclusion  of 
these  studies  is,  that  them  are  certain  application  areas  where  an  ANN  implemented  as  an 
analog  VLSI  chip  could  perform  critical  functions  as  a  component  of  a  hybrid  system.  The 
advantage  of  such  an  approach  is  the  high  speed,  low  weight  and  power  consumption,  and 
’’built-in”  fault-tolerance  of  the  ANN,  which  is  especially  attractive  for  long-term  missions, 
where  component  failures  have  to  be  expected  but  a  repair  is  not  possible. 

In  order  to  extend  the  scope  of  ANN  models  and  application  areas  considered  so  far, 
a  new  project  in  cooperation  with  D.  Palumbo  and  C.  Jorgensen  (NASA  Ames  Research 
Center)  investigates  new  unsupervised  learning  methods  for  nonlinear,  adaptive  control  with 
potential  use  in  autonomous  systems.  From  a  fault-tolerant  systems  perspective,  there  are 
three  different  issues  that  can  be  addressed,  a  failure  of  the  controller,  a  failure  of  the  plant, 
and  a  sudden  change  in  the  environmental  conditions  or  external  perturbations.  For  example, 
an  adaptive  ANN  controlling  a  robot  arm  should  be  able  to  recover  from  internal  faults,  from 
a  sudden  increase  in  the  friction  of  an  arm  joint,  or  from  a  new  obstacle  in  the  environment 
by  updating  and  adjusting  its  control  function.  The  main  problem  is  to  maintain  a  degree 
of  plasticity  of  the  ANN,  but  still  to  guarantee  an  overall  stable  behavior. 

Work  continues  on  another  project  in  collaboration  with  C.  Jeffries  (Clemson  University) 
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that  investigates  a  new  associative  memory  architecture  with  high  order  feedback  connec¬ 
tions.  One  application  that  has  been  successfully  implemented  is  the  error  correction  of  block 
codes  transmitted  over  a  noisy  channel  in  a  digital  communication  system.  Simulations  have 
shown  that  the  use  of  the  associative  memory  reduces  the  achievable  bit  error  rate  by  up 
to  two  orders  of  magnitude  in  comparison  to  conventional  methods.  We  are  currently  work¬ 
ing  on  a  general  stability  proof  of  the  systems  which  would  extend  earlier  results  that  exist 
only  for  special  cases  of  the  ’’neuron”  transfer  function.  Furthermore,  we  plan  to  investigate 
extensions  of  the  architecture  and  different  signal  processing  applications. 


Joel  Saltz 

The  activities  described  below  were  done  in  collaboration  with  S.  Berryman,  K.  Crowley, 
R.  Das,  S.  Gupta  (Old  Dominion  University),  S.  Hiranandani,  D.  Mavriplis,  S.  Petiton 
(ETCA,  France),  R.  Ponnusamy,  J.  Scroggs,  and  J.  Wu  (Yale  University). 

We  have  developed  a  prototype  compiler  which  takes  as  input  a  Fortran  77  program 
enhanced  with  specifications  for  distributing  data.  The  compiler  outputs  a  message  passing 
program  that  runs  on  a  distributed  memory  computer.  The  runtime  support  for  this  compiler 
is  the  PARTI  library  of  primitives  designed  to  efficiently  support  irregular  patterns  of  dis¬ 
tributed  array  accesses  and  irregular  distributed  array  partitions.  This  compiler  was  tested 
on  several  NASA  kernels,  and  the  performance  of  the  resulting  codes  was  benchmarked  on 
the  iPSC/860.  A  description  of  the  compiler  and  the  associated  computational  experiments 
may  be  found  in  ICASE  Reports  90-59  and  91-13. 

We  have  developed  ana  distributed  two  new  3ets  of  runtime  primitives  designed  to  help 
users  handle  sparse  and  unstructured  computations.  The  new  primitives  incorporate  a  num¬ 
ber  of  new  insights  we  have  had  about  sparse  and  unstructured  computations.  Our  new 
primitives  carry  out  preprocessing  that  makes  it  straightforward  to  produce  parallelized 
loops  that  are  virtually  identical  in  form  to  the  original  sequential  loops.  The  importance  of 
this  is  that  it  will  be  possible  to  generate  the  same  quality  object  code  on  the  nodes  of  the 
distributed  memory  machine  as  could  be  produced  by  the  sequential  program  running  on  a 
single  node. 

Our  new  primitives  incorporate  hash  tables;  hash  tables  are  used  to  allow  us  to  recognize 
situations  in  which  a  single  off-processor  distributed  array  reference  is  used  several  times. 
The  new  primitives  only  fetch  a  single  copy  of  each  unique  off-processor  distributed  array 
reference.  The  new  primitives  also  allow  us  to  make  much  finer  distinctions  concerning  which 
new  data  need  to  be  obtained  from  other  processors.  In  our  system,  schedules  are  generated 
to  specify  which  distributed  array  elements  are  to  be  moved  between  processors.  In  the 
older  versions  of  the  primitives,  we  had  to  generate  each  schedule  from  scratch;  we  can  now 
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produce  a  schedule  that  allows  us  to  obtain  only  those  off-processor  elements  that  have  not 
been  previously  encountered  by  a  specified  set  of  other  schedules.  Some  of  the  methods  used 
to  accomplish  this  are  described  in  ICASE  Report  90-33  and  a  full  description  of  this  work 
will  be  forthcoming. 

The  new  set  of  primitives  have  been  used  to  port  a  3-D  unstructured  mesh  Euler  solver. 
This  work  has  spurred  many  improvements  in  the  optimizations  carried  out  by  our  primitives; 
we  have  seen  a  reduction  in  communication  time  in  a  problem  solved  on  a  50,000  node  grid, 
from  288  seconds  to  82  seconds  (compared  to  a  computation  time  of  151  seconds).  In  this 
code,  the  total  cost  of  all  preprocessing  was  less  than  2  seconds. 

Since  the  form  of  the  sequential  code  and  the  parallelized  code  was  virtually  identical,  we 
did  not  expect  the  parallelization  process  to  introduce  any  new  inefficiencies  beyond  those 
exacted  by  the  preprocessing  and  by  the  calls  to  the  primitives.  On  a  smaller  problem,  we 
compared  the  parallel  code  running  on  a  single  node  with  the  sequential  code  and  found  only 
a  2  percent  performance  degradation. 

We  have  also  made  progress  in  the  area  of  primitives  for  structured  mesh  adaptive  codes 
and  for  block  structured  codes.  In  block  structured  codes  (such  as  CFL3D  and  GASP), 
solvers  are  carried  out  on  a  varying  number  of  coupled  Cartesian  meshes.  The  meshes  may 
be  different  sizes,  and  each  mesh  needs  to  be  mapped  to  a  subset  of  a  machine’s  processors. 

In  structured  adaptive  mesh  problems,  the  relationship  between  the  meshes  that  define 
a  problem  change  dynamically.  In  many  of  these  codes,  each  computational  phase  involves 
only  a  subset  of  the  declared  distributed  arrays.  Since  we  want  to  achieve  reasonable  load 
balance  in  each  computational  phase,  we  need  to  be  able  to  partition  the  various  distributed 
arrays  in  an  overlapping  manner.  The  distributed  array  mappings  are  determined  at  runtime 
and  may  be  dynamic. 

The  primitives  aimed  at  block  structured  and  at  structured  adaptive  problems  incorporate 
geometrical  information.  These  primitives  allow  users  to  dynamically  initialize  multidimen¬ 
sional  distributed  array  mappings  as  well  as  to  transfer  data  between  irregularly  distributed 
multidimensional  arrays.  We  are  currently  working  on  a  version  of  the  primitives  that  will 
support  dynamic  array  remapping.  Some  of  this  work  is  described  in  ICASE  Report  90-41. 

The  PARTI  runtime  support  mechanisms  are  also  playing  an  influential  role  in  a  variety 
of  distributed  memory  compiler  development  projects  in  the  United  States  and  abroad.  The 
primitives  have  been  distributed  to  the  compiler  groups  at  Rice  University,  University  of 
Colorado,  Cornell  University,  University  of  Washington,  and  University  of  Vienna  among 
others. 

We  are  also  extending  this  work  to  SIMD  multiprocessors.  We  are  developing  a  set  of 
MasPar  PARTI  primitives  for  the  MP-1  built  by  MasPar  Inc.  and  intend  to  use  them  to  port 
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the  3-D  unstructured  Euler  cude  discussed  earlier.  We  have  also  carried  out  a  comprehensive 
benchmarking  project  involving  a  comparison  of  the  effects  of  irregular  communications  and 
of  problem  mapping  on  the  CM-2  built  by  Thinking  Machines  Inc.  This  study  was  reported 
in  ICASE  Report  91-12. 


Joel  Saltz  and  D.  Mavriplis 

The  thin  layer  Navier-Stokes  equations  are  solved  for  two-dimensional  airfoil  problems  by 
preconditioned  conjugate  gradient-like  iterative  methods.  In  past  work,  V.  Venkatakrishnan 
(Computer  Science  Corporation,  Moffett  Field,  CA)  has  found  the  GMRES  with  incomplete 
LU  factorization  as  a  preconditioner  is  an  excellent  scheme  for  solving  linear  equations  arising 
from  two-dimensional  airfoil  flow  calculations  on  structured  meshes.  We  have  shown  that 
these  methods  are  equally  applicable  to  unstructured  grids. 

In  collaboration  with  Venkatakrishnan,  we  have  investigated  a  variety  of  issues  that  arise 
when  mapping  such  iterative  methods  to  vector  parallel  machines  and  have  characterized 
their  performance.  The  inner  loop  of  our  solver  involves  computation  of  upper  and  lower 
sparse  triangular  systems.  These  triangular  matrix  solutions  can  constitute  an  appreciable 
percentage  of  the  operations  required  in  the  iterative  portion  of  preconditioned  Krylov  space 
algorithms.  It  is  consequently  essential  to  efficiently  compute  these  sparse  triangular  solves. 
We  have  performed  the  optimizations  required  to  vectorize  and  parallelize  these  operations 
on  an  eight  processor  Cray  Y/MP. 

While  we  achieved  satisfactory  results  in  vectorizing  the  sparse  triangular  solves,  our 
vector  parallel  results  were  still  somewhat  disappointing.  Another  approach  is  to  partition 
a  domain  into  P  subdomains  and  carry  out  an  ILU  factorization  over  each  sub-domain. 
The  convergence  properties  of  such  a  domain  decomposition  depends  on  precisely  how  the 
domain  is  subdivided.  When  we  partition  the  domain  in  one  of  our  model  problems  in  a 
judicious  manner,  we  are  able  to  obtain  speedups  of  approximately  4.4  compared  to  a  single 
domain,  vectorized  code.  There  is  an  overhead  of  approximately  20%  due  to  extra  operations 
introduced;  most  of  the  rest  of  the  overhead  appears  to  be  due  to  the  fact  that  partitioning 
a  problem  into  subdomains  reduces  vector  lengths.  We  have  investigated  several  different 
methods  for  partitioning  unstructured  meshes  to  obtain  good  convergence  and  good  load 
balance,  and  to  maximize  vector  lengths.  We  have  also  attempted  to  model  the  Cray  Y /MP 
performance  so  that  we  can  extend  our  results  to  deal  with  future  more  highly  parallel  vector 
machines.  These  results  will  appear  in  the  proceedings  of  the  SIAM  Conference  on  Parallel 
Processing  and  as  an  JCASE  report. 
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Sutanu  Sarkar 

We  are  engaged  in  the  direct  simulation  and  Reynolds  stress  modeling  of  turbulent  flows. 
The  case  of  compressible  homogeneous  shear  flow  has  been  simulated  in  collaboration  with  G. 
Erlebacher  and  M.  Y.  Hussaini.  The  simulations  reveal  a  number  of  features  associated  with 
compressibility,  for  example,  a  reduction  in  the  growth  rate  of  the  turbulent  kinetic  energy, 
the  augmentation  of  the  turbulent  dissipation  rate,  and  a  quasi-equilibrium  structure  of  the 
dilatational  component  of  the  velocity  field.  We  have  shown  clearly  that  the  reduced  growth 
of  turbulent  kinetic  energy  in  the  compressible  case  is  associated  with  the  compressible 
dissipation  and  the  pressure-dilatation.  We  have  also  obtained  results  regarding  probability 
density  functions  and  higher-order  moments  from  the  simulations.  These  results  indicate  that 
local  fluid  compressions  are  more  likely  than  local  expansions,  and  that  the  compressible 
velocity  component  is  more  intermittent  than  the  vortical  velocity  component.  Reynolds 
stress  modeling  for  high-speed  flows  is  continuing  in  collaboration  with  C.  Speziale  and  T. 
Gatski  (Fluid  Mechanics  Division,  LaRC).  We  now  have  a  turbulence  closure  for  the  averaged 
Navier-Stokes  equations,  the  Reynolds  stress  equations  and  the  turbulent  dissipation  rate. 
The  new  compressibility  models  and  the  SSG  pressure-strain  model  have  been  incorporated 
into  this  closure.  Application  of  the  Reynolds  stress  closure  to  the  supersonic  shear  layer 
and  the  boundary  layer  have  produced  promising  results. 


Sharon  O.  Seddougui 

Work  is  continuing  on  the  spatial  inviscid  instability  problem  for  Gortler  vortices  in 
a  compressible  fluid.  The  governing  equations  in  the  limit  of  large  Gortler  number  were 
solved  numerically  in  conjunction  with  Andrew  Dando  (University  of  Manchester)  to  find 
the  growth  rate  of  the  imposed  disturbance.  We  see  the  appearance  of  two  distinct  modes  of 
instability,  one  of  which  is  almost  continuous  for  large  wavenumber.  Asymptotic  solutions 
have  been  obtained  for  large  Mach  number  and/or  large  wavenumber.  These  are  used  to 
explain  the  near-linking  of  the  modes  which  occur.  A  joint  report  with  Andrew  Dando  is 
being  prepared. 

The  extension  of  the  work  by  Cowley  and  Hall  (ICASE  Report  No.  88-72)  on  the 
instability  of  hypersonic  flow  past  a  wedge  to  include  the  effects  of  nonlinearity  has  been 
completed  jointly  with  Andrew  P.  Bassom  (University  of  Exeter).  We  are  in  the  process 
of  writing  a  report  on  the  results.  We  found  that  the  disturbance  is  supercritical  i.e.  the 
nonlinear  effects  are  stabilizing  and  that  a  stable  nonzero  equilibrium  amplitude  is  possible. 

The  work  of  MacKerrell  (Ph.D.  thesis,  University  of  Exeter),  which  investigates  the  effect 
of  compliant  walls  on  the  hydrodynamic  stability  of  a  boundary  layer  flow,  has  been  continued 
to  investigate  the  effects  of  nonlinearity.  The  linear  stability  analysis  shows  that  the  growth 
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rate  of  a  disturbance  may  be  reduced  by  having  a  flexible  boundary.  It  is  found  that  the 
nonlinear  effects  are  stabilizing  and  that  there  is  the  possibility  of  an  equilibrium  solution  if 
the  flow  is  linearly  unstable.  These  results  are  being  prepared  for  an  ICASE  Report. 


Charles  G.  Speziale 

A  compressible  second-order  closure  model  for  supersonic  turbulent  flows  has  been  devel¬ 
oped  in  collaboration  with  S.  Sarkar.  The  model  has  been  tested  in  a  variety  of  homogeneous 
and  inhomogeneous  supersonic  turbulent  flows  and  comparisons  with  the  results  from  physi¬ 
cal  and  numerical  experiments  are  quite  favorable.  Further  tests,  and  possible  refinements  of 
the  model,  are  planned  for  future  research.  Additional  tests  of  a  recently  developed  near-wall 
two-equation  turbulence  model  (ICASE  Report  No.  90-46)  were  conducted  in  collaboration 
with  S.  Thangam  and  R.  Abid  (Vigyan  Research  Associates).  The  model  was  shown  to  per¬ 
form  extremely  well  in  two  independent  near- wall  turbulent  flows  involving  adverse  pressure 
gradients  and  flow  separation.  A  computational  study  of  the  performance  of  two-equation 
turbulence  models  for  turbulent  flow  past  a  backward-facing  step  was  also  conducted  in 
collaboration  with  S.  Thangam.  It  was  found  that  previously  reported  inaccuracies  in  the 
predictions  of  the  standard  K  —  e  model  were  overestimated  due  to  inadequate  numerical 
resolution  of  the  flow  field.  With  proper  numerical  resolution  -  and  with  the  use  of  an 
anisotropic  eddy  viscosity  -  the  K  —  e  model  was  shown  to  yield  excellent  predictions  for 
the  backstep  problem. 

Independent  work  has  continued  with  T.  Gatski  (Fluid  Mechanics  Division,  LaRC)  and 
P.  Durbin  (CTR,  Stanford  University)  on  modeling  the  tensor  dissipation  rate  of  turbulence. 
The  work  with  P.  A.  Durbin  has  cast  grave  suspicions  on  the  general  validity  of  Kolmogorov’s 
hypothesis  of  local  isotropy.  This  hypothesis  -  which  has  served  as  a  cornerstone  for  most 
turbulence  models  -  postulates  that,  at  sufficiently  high  Reynolds  numbers,  the  small  scales 
of  an  anisotropic,  inhomogeneous  turbulence  are  approximately  isotropic.  If  this  hypothesis 
is  not  valid,  models  for  the  tensor  dissipation  rate  could  play  a  crucial  role  in  the  development 
of  improved  turbulence  models.  Consequently,  future  research  is  planned  in  this  direction. 


Eitan  Tadmor 

We  continue  our  study  of  the  isentropic  gasdynamics  equations  using  ’appropriate’  kinetic 
models,  which  enable  -  by  averaging  compactness  arguments  -  derivation  of  new  results  on 
macroscopic  regularity.  Also,  a  regularized  Chapman-  Enskog  expansion  for  scalar  conserva¬ 
tion  laws  is  studied  in  ICASE  Report  No.  90-68,  where  a  bounded  (relaxation)  approximation 
of  the  linearized  collision  operator  is  used  for  high  wave  numbers. 


22 


We  continue  the  development  of  a  convergence  analysis  for  the  approximate  solution  to 
nonlinear  hyperbolic  equations. 

We  are  completing  our  study  of  entropy  stable  difference  approximations  to  nonlinear 
systems  of  conservation  laws. 

Finally,  we  began  to  study  various  edge  detectors  based  on  the  (truncated)  Hilbert  trans¬ 
form,  in  the  context  of  spectral  data. 


Saleh  Tanveer 

Progress  has  been  made  in  several  areas  of  fluid  dynamics  and  crystal  growth  in  the 
presence  of  convection. 

In  the  area  of  water  waves  and  the  classical  Rayleigh-Taylor  flow,  we  found  that  the 
characterization  of  the  flow  field  in  terms  of  singularities  in  the  unphysical  domain  was  useful 
from  a  computational  viewpoint.  Such  an  approach  had  the  potential  of  calculating  highly 
convoluted  interface  and  explaining  the  interaction  of  bubbles  and  spikes  for  the  Rayleigh- 
Taylor  flow  in  terms  of  interaction  of  singularities.  This  work  is  described  in  ICASE  Report 
91-21. 

In  the  initial  value  problem  for  an  interface  in  a  Hele-Shaw  cell,  a  problem  that  is  ill  posed 
in  the  absence  of  surface  tension,  we  found  that  a  convenient  first  step  in  the  understanding  of 
the  regularizing  effect  of  small  nonzero  surface  tension  was  to  embed  the  zero  surface  tension 
problem  into  a  well  posed  problem  by  extending  the  spatial  domain  to  the  unphysical  region. 
Small  surface  tension  effects  can  then  be  studied  perturbatively.  We  found  this  approach 
to  be  fruitful  in  explaining  highly  convoluted  interfacial  features  that  are  very  sensitive  to 
initial  conditions.  Our  approach  may  be  general  with  possible  application  to  Euler  flows 
and  understanding  of  effect  of  small  viscosity.  An  ICASE  report  describing  the  Hele-Shaw 
dynamics  is  currently  being  prepared. 

For  the  problem  of  crystal  growth  in  a  cylindrical  geometry  with  a  forced  horizontal  ther¬ 
mal  gradient,  we  have  some  regions  of  the  fifteen  parameter  space  with  different  dynamics. 
This  has  all  been  done  for  small  values  of  heat  transfer  through  the  side  walls.  We  have 
scaling  results  on  the  relative  size  of  convection  terms  near  the  free  interface  for  different 
ranges  of  solutal  and  thermal  Rayleigh  numbers.  Further  work  is  in  progress  to  uncover  all 
other  interesting  limits. 


Siva  Thangam 

Computational  and  experimental  investigations  of  turbulent  separated  flow  past  a  rear¬ 
ward  facing  step  are  underway  to  examine  the  efficacy  of  turbulence  models.  The  first  phase 
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of  this  collaborative  effort  involving  C.  Speziale,  M.Y.  Hussaini,  and  S.  Kjelgaard  (Fluid 
Mechanics  Division,  LaRC)  has  primarily  been  concerned  with  the  use  of  Reynolds  stress 
closure  models.  The  full  Navier-Stokcs  equations  and  the  transport  equations  for  the  tur¬ 
bulence  kinetic  energy  and  dissipation  are  solved  using  different  wall  boundary  conditions 
and  two-equation  turbulence  models  with  isotropic  and  anisotropic  eddy  viscosity  represen¬ 
tations.  The  computations  are  performed  using  a  staggered  mesh,  finite-volume  algorithm. 
The  results  indicate  that  the  characteristics  of  the  flow  field  are  sensitive  to  the  type  of 
turbulence  model,  the  representation  of  the  wall  conditions,  and  the  inflow  conditions.  It  is 
demonstrated  that  if  the  computations  were  performed  such  that  the  inflow  and  the  bound¬ 
ary  conditions  essentially  replicate  those  in  the  experiments  very  good  agreement  between 
the  experimental  and  computational  results  could  be  achieved.  During  this  phase,  the  ex¬ 
perimental  setup  was  also  redesigned  with  an  improved  inlet  channel  and  a  substantially 
longer  downstream  test  section.  During  the  second  phase,  experiments  will  be  performed 
for  a  step  to  channel  height  ratio  of  1:2  at  several  inflow  conditions.  Computations  also 
will  be  performed  simultaneously  using  two-equation  turbulence  models  with  isotropic  and 
anisotropic  eddy  viscosity  representations  and  the  predictions  will  be  compared  with  those 
from  the  experiments.  In  addition,  large-eddy  simulation  of  the  turbulent  flow  past  the  1:2 
rearward  facing  step  will  be  performed  to  evaluate  the  performance  of  the  Reynolds  stress 
closure  models.  The  issues  involving  the  appropriate  representation  for  the  subgrid  scale 
stresses  will  be  addressed. 


John  Van  Rosendale 

Work  done  in  collaboration  with  Andrea  Ove:  man  (Analysis  and  Computation  Division, 
LaRC)  found  that  spectral  methods  require  frequent  global  communication,  raising  a  fun¬ 
damental  barrier  to  their  use  on  parallel  architectures.  Global  communication  is  needed  in 
both  the  spectral  derivative  computation,  and  in  the  implicit  solves.  One  can  do  the  implicit 
solves  in  several  different  ways,  but  one  of  the  best  alternatives  is  to  use  finite  element  pre¬ 
conditioning,  as  pioneered  by  Deville  and  Mund.  The  preconditioning  finite  element  system 
can  in  turn  be  solved  by  a  multigrid  algorithm.  A  well  tuned  multigrid  algorithm  adequately 
solves  the  preconditioning  finite  element  system  in  one  V-cycle. 

To  explore  the  issue  of  communication  in  spectral  methods,  we  implemented  a  three 
dimensional  implicit  spectral  code  for  the  variable  coefficient  Helmholtz  equation  on  the 
Intel  iPSC/860.  This  equation  arises  frequently  in  acoustics  and  in  the  Uzawa  formulation 
of  the  incompressible  Navier  Stokes  equation.  Our  Fourier-spectral  code,  based  on  a  V- 
cycle  point-relaxation  mu'ltigrid  algorithm  achieves  utilization  of  50%  on  the  ICASE  32  node 
iPSC/860,  for  a  64  x  64  x  64  grid;  finer  grids  yield  higher  utilizations. 
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Chebyshev-spectral  grids  are  more  problematic,  since  either  plane-relaxation  based  multi¬ 
grid,  or  the  new  robust  point-relaxation  algorithm  is  required.  We  implemented  a  semicoars¬ 
ening  plane-relaxation  algorithm,  using  concurrent  relaxation  between  planes  to  increase  par¬ 
allelism.  This  approach  leads  to  perfect  load  balance,  and  is  well  suited  to  machines  like  the 
iPSC/860,  having  a  modest  number  of  processors.  We  achieve  typical  utilizations  of  about 
65%.  However,  this  multigrid  algorithm  solves  the  preconditioning  finite  element  systems 
more  thoroughly  than  necessary,  so  the  new  point-relaxation  based  multigrid  algorithms  may 
be  preferable.  Exploration  of  this  alternative  is  currently  in  progress. 
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REPORTS  AND  ABSTRACTS 


Bernard,  Peter  S.,  and  Charles  G.  Speziale:  Bounded  energy  states  in  homogeneous  turbulent 
shear  flow  -  An  alternative  view.  ICASE  Report  No.  90-66,  October  2,  1990,  35  pages.  To 
be  submitted  to  Physics  of  Fluids  A. 

The  equilibrium  structure  of  homogeneous  turbulent  shear  flow  is  investigated  from  a 
theoretical  standpoint.  Existing  turbulence  models,  in  apparent  agreement  with  physical 
and  numerical  experiments,  predict  an  unboundedexponential  time  growth  of  the  turbulent 
kinetic  energy  and  dissipation  rate;  only  the  anisotropy  tensor  and  turbulent  time  scale 
reach  a  structural  equilibrium.  It  is  shown  that  if  vortex  stretching  is  accounted  for  in  the 
dissipation  rate  transport  equation,  then  there  can  exist  equilibrium  solutions, with  bounded 
energy  states,  where  the  turbulence  production  is  balanced  by  its  dissipation.  Illustrative 
calculations  are  presented  for  a  k  —  e  model  modified  to  account  for  vortex  stretching.  The 
calculations  indicate  an  initial  exponential  time  growth  of  the  turbulent  kinetic  energy  and 
dissipation  rate  for  elapsed  times  that  are  as  large  as  those  considered  in  any  of  the  previously 
conducted  physical  or  numerical  experiments  on  homogeneous  shear  flow.  However, vortex 
stretching  eventually  takes  over  and  forces  a  production-equals-dissipation  equilibrium  with 
bounded  energy  states.  The  validity  of  this  result  is  further  supported  by  an  independent 
theoretical  argument.  It  is  concluded  that  the  generally  accepted  structural  equilibrium  for 
homogeneous  shear  flow  with  unbounded  component  energies  is  in  need  of  re-examination. 


Nicol,  David  M.,  and  Scott  E.  Riffe:  A  “ conservative ”  approach  to  parallelizing  the  sharks 
world  simulation.  ICASE  Report  No.  90-67,  October  2,  1990,  14  pages.  To  appear  1990 
Winter  Simulation  Conf.  Proceedings. 


This  paper  describes  how  we  parallelized  a  benchmark  problem  for  parallel  simulation, 
the  Sharks  World.  The  solution  we  describe  is  conservative,  in  the  sense  that  no  state 
information  in  saved,  and  no  “rollbacks”  occur.  Our  approach  illustrates  both  the  principal 
advantage  and  principal  disadvantage  of  conservative  parallel  simulation.  The  advantage 
is  that  by  exploiting  lookahead  we  find  an  approach  that  dramatically  improves  the  serial 
execution  time,  and  ako  achieves  excellent  speedups.  The  disadvantage  is  that  if  the  model 
rules  are  changed  in  such  a  way  that  the  lookahead  is  destroyed,  it  is  difficult  to  modify  the 
solution  to  accommodate  the  changes. 


Schochet,  Steven,  and  Eitan  Tadmor:  Regularized  Chapman-Enskog  expansion  for  scalar 
conservation  laws.  ICASE  Report  No.  90-68,  October  2,  1990,  19  pages.  Submitted  to  the 
Archive  for  Mechanics  and  Analysis. 

Rosenau  [Phys.  Rev.  A,  40  (1989),  pp.  7193-6]  has  recently  proposed  a  regularized 
version  of  the  Chapman-Enskog  expansion  of  hydrodynamics.  This  regularized  expansion 
resembles  the  usual  Navier-Stokes  viscosity  terms  at  law  wavenumbers,  but  unlike  the  latter, 
it  has  the  advantage  of  being  a  bounded  macroscopic  approximation  to  the  linearized  collision 
operator. 

In  this  paper  we  study  the  behavior  of  Rosenau  regularization  of  the  Chapman-Enskog 
expansion  (R-C-E)  in  the  context  of  scalar  conservation  laws.  We  show  that  this  R-C-E 
model  retains  the  essential  properties  of  the  usual  viscosity  approximation,  e.g.,  existence 
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of  travelling  waves,  monotonicity,  upper- Lipschitz  continuity,  etc.,  and  at  the  same  time,  it 
sharpens  the  standard  viscous  shock  layers.  We  prove  that  the  regularized  R-C-E  approx¬ 
imation  converges  to  the  underlying  inviscid  entropy  solution  as  its  mean-free  path  £  |  0, 
and  we  estimate  the  convergence  rate. 


Mehrotra,  Piyush,  and  John  Van  Rosendale:  Programming  distributed  memory  architectures 
using  Kali.  ICASE  Report  No.  90-69,  October  3,  1990,  23  pages.  Proceedings  of  the  Work¬ 
shop  on  Languages  and  Compilers  for  Parallel  Computers,  MIT/Pitman  Press,  Fall  1990. 


Programming  nonshared  memory  systems  is  more  difficult  than  programming  shared 
memory  systems,  in  part  because  of  the  relatively  low  level  of  current  programming  envi¬ 
ronments  for  such  machines.  This  paper  presents  a  new  programming  environment,  Kali, 
which  provides  a  global  name  space  and  allows  direct  access  to  remote  data  values.  In  order 
to  retain  efficiency,  Kali  provides  a  system  of  annotations,  allowing  the  user  to  control  those 
aspects  of  the  program  critical  to  performance,  such  as  data  distribution  and  load  balanc¬ 
ing.  This  paper  describes  the  primitives  and  constructs  provided  by  our  language,  and  also 
discusses  some  of  the  issues  raised  in  translating  a  Kali  program  for  execution  on  distributed 
memory  systems. 


Koelbel,  Charles,  and  Piyush  Mehrotra:  Compiling  global  name-space  programs  for  dis¬ 
tributed  execution.  ICASE  Report  No.  90-70,  October  4,  1990,  36  pages.  Submitted  to 
IEEE  Transactions  on  Parallel  and  Distributed  Computing,  Special  Issue  June  1991. 


Distributed  memory  machines  do  not  provide  hardware  support  for  a  global  address 
space.  Thus  programmers  are  forced  to  partition  the  data  across  the  memories  of  the  ar¬ 
chitecture  and  use  explicit  message  passing  to  communicate  data  between  processors.  In 
this  paper  we  focus  on  the  compiler  support  required  to  allow  programmers  to  express  their 
algorithms  using  a  global  name-space.  We  present  a  general  method  for  analysis  of  a  high 
level  source  program  and  its  translation  to  a  set  of  independently  executing  tasks  communi¬ 
cating  via  messages.  If  the  compiler  has  enough  information,  this  translation  can  be  carried 
out  at  compile-time.  Otherwise  run-time  code  is  generated  to  implement  the  required  data 
movement.  We  describe  the  analysis  required  in  both  situations  and  present  the  performance 
of  the  generated  code  on  the  Intel  iPSC/2. 


Hall,  Philip,  and  Nicola  J.  Horseman:  The  inviscid  secondary  instability  of  fully  nonlinear 
longitudinal  vortex  structures  in  growing  boundary  layers.  ICASE  Report  No.  90-71,  Octo¬ 
ber  10,  1990,  41  pages.  Submitted  to  Journal  of  Fluid  Mechanics. 


The  inviscid  instabilty  of  a  longitudinal  vortex  structure  within  a  steady  boundary  layer  is 
investigated.  The  instability  has  wavelength  comparable  with  the  boundary  layer  thickness  so 
that  a  quasi-parallel  approach  to  the  instability  problem  can  be  justified.  The  generalization 
of  the  Rayleigh  equation  to  such  a  flow  is  obtained  and  solved  for  the  case  when  the  vortex 
structure  is  induced  by  curvature.  Two  distinct  modes  of  instability  are  found;  these  modes 
correspond  with  experimental  observations  on  the  breakdown  process  for  Gortler  vortices. 
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Bassom,  Andrew  P.,  and  Philip  Hall:  Vortex  instabilities  in  SD  boundary  layers:  The  re¬ 
lationship  between  Gortler  and  crossflow  vortices.  ICASE  Report  No.  90-72,  October  11, 
1990,  76  pages.  Submitted  to  Journal  of  Fluid  Mechanics. 


The  inviscid  and  viscous  stability  problems  are  addressed  for  a  boundary  layer  which  can 
support  both  Gortler  and  Crossflow  vortices.  The  change  in  structure  of  Gortler  vortices 
is  found  when  the  parameter  representing  the  degree  of  three-dimensionality  of  the  basic 
boundary  layer  flow  under  consideration  is  increased.  It  is  shown  that  Crossflow  vortices 
emerge  naturally  as  this  parameter  is  increased  and  ultimately  become  the  only  possible 
vortex  instability  of  the  flow.  It  is  shown  conclusively  that  at  sufficiently  large  values  of 
the  crossflow  there  are  no  unstable  Gortler  vortices  present  in  a  boundary  layer  which,  in 
the  zero  crossflow  case,  is  centrifugally  unstable.  The  results  suggest  that  in  many  practical 
applications  Gortler  vortices  cannot  be  a  cause  of  transition  because  they  are  destroyed  by 
the  3-D  nature  of  the  basic  state.  In  swept  wing  flows  the  Gortler  mechanism  is  probably  not 
present  for  typical  angles  of  sweep  of  about  20  degrees.  Some  discussion  of  the  receptivity 
problem  for  vortex  instabilities  in  weakly  3-D  boundary  layers  is  given;  it  is  shown  that 
inviscid  modes  have  a  coupling  coefficient  marginally  smaller  than  those  of  the  fastest  growing 
viscous  modes  discussed  recently  by  Denier,  Hall  and  Seddougui  (1990).  However  the  fact 
that  the  growth  rates  of  the  inviscid  modes  are  the  largest  in  most  situations  means  that 
they  are  probably  the  most  likely  source  of  transition. 


Abarbanel,  Saul,  and  David  Gottlieb:  Spurious  frequencies  as  a  result  of  numerical  bound¬ 
ary  treatments.  ICASE  Report  No.  90-73,  October  12,  1990,  28  pages.  To  appear  in  the 
Proceedings  of  the  Third  International  Conference  on  Hyperbolic  Systems  (invited  paper). 


The  stability  theory  for  finite  difference  Initial  Boundary- Value  approximations  to  sys¬ 
tems  of  hyperbolic  partial  differential  equations  states  that  the  exclusion  of  eigenvalues  and 
generalized  eigenvalues  is  a  sufficient  condition  for  stability.  The  theory,  however,  does  not 
discuss  the  nature  of  numerical  approximations  in  the  presence  of  such  eigenvalues. 

In  fact,  as  was  shown  previously  [1],  for  the  problem  of  vortex  shedding  by  a  2-D  cylinder 
in  subsonic  flow,  stating  boundary  conditions  in  terms  of  the  primitive  (non-characteristic) 
variables  may  lead  to  such  eigenvalues,  causing  perturbations  that  decay  slowly  in  space  and 
remain  periodic  time.  Characteristic  formulation  of  the  boundary  conditions  avoided  this 
problem. 

In  this  paper,  we  report  on  a  more  systematic  study  of  the  behavior  of  the  (linearized) 
one-dimensional  gas  dynamic  equations  under  various  sets  of  oscillation-inducing  “legal” 
boundary  conditions. 

Gottlieb,  D.,  and  C.  L.  Streett:  Quadrature  imposition  of  compatability  conditions  in  Cheby- 
shev  methods.  ICASE  Report  No.  90-74,  October  12,  1990,  18  pages.  Submitted  to  Journal 
of  Scientific  Computing. 


Often,  in  solving  an  elliptic  equation  with  Neumann  boundary  conditions,  a  compatibility 
condition  has  to  be  imposed  for  well-posedness.  Thiscondition  involves  integrals  of  the  forcing 
function. 

When  pseudospectral  Chebyshev  methods  are  used  to  discretize  the  partial  differential 
equation,  these  integrals  have  to  be  approximated  by  an  appropriate  quadrature  formula. 
The  Gauss-Chebyshev  (or  any  variant  of  it,  like  the  Gauss-Lobatto)  formula  can  not  be  used 
here  since  the  integrals  under  consideration  do  not  include  the  weight  function.  A  natural 
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candidate  to  be  used  in  approximating  the  integrals  is  the  Clenshaw- Curtis  formula,  however 
we  show  in  this  paper  that  this  is  the  wrong  choice  and  it  may  lead  to  divergence  if  time 
dependent  methods  are  used  to  march  the  solution  to  ?  ueady  state. 

We  develop,  in  this  paper,  the  correct  quadrature  formula  for  these  problems.  This 
formula  takes  into  account  the  degree  of  the  polynomials  involved.  We  show  that  this 
formula  leads  to  a  well  conditioned  Chebyshev  approximation  to  the  differential  equations 
and  that  the  compatibility  condition  is  automatically  satisfied. 


Roe,  P.  L.:  “Optimum”  upwind  advection  on  a  triangular  mesh.  ICASE  Report  No.  90-75, 
October  15,  22  pages.  To  be  submitted  to  Journal  of  Computational  Physics. 


For  advection  schemes  based  on  fluctuation  splitting,  a  design  criterion  of  optimising  the 
time  step  leads  to  linear  schemes  that  coincide  with  those  designed  for  least  truncation  error. 
A  further  stage  of  optimising  the  time  step  using  a  non-linear  positivity  criterion,  leads  to 
considerable  further  gains  in  resolution. 


Erlebacher,  G.,  M.  Y.  Hussaini,  C.  G.  Speziale,  and  T.  A.  Zang:  Toward  the  large-eddy 
simulation  of  compressible  turbulent  flows.  ICASE  Report  No.  90-76,  October  16,  1990,  45 
pages.  Submitted  to  the  Journal  of  Fluid  Mechanics. 


New  subgrid-scale  models  for  the  large-eddy  simulation  of  compressible  turbulent  flows 
are  developed  and  tested  based  on  the  Favre-filtered  equations  of  motion  for  an  ideal  gas.  A 
compressible  generalization  of  the  linear  combination  of  the  Smagorinsky  model  and  scale- 
similarity  model,  in  terms  of  Favre-filtered  fields,  is  obtained  for  the  subgrid-scale  stress 
tensor.  An  analogous  thermal  linear  combination  model  is  also  developed  for  the  subgrid- 
scale  heat  flux  vector.  The  two  dimensionless  constants  associated  with  these  subgrid-scale 
models  are  obtained  by  correlating  with  the  results  of  direct  numerical  simulations  of  com¬ 
pressible  isotropic  turbulence  performed  on  a  963  grid  using  Fourier  collocation  methods. 
Extensive  comparisons  between  the  direct  and  modeled  subgrid-scale  fields  are  provided  in 
order  to  validate  the  models.  A  large-eddy  simulation  of  the  decay  of  compressible  isotropic 
turbulence  -  conducted  on  a  coarse  323  grid  -  is  shown  to  yield  results  that  are  in  excellent 
agreement  with  the  fine  grid  direct  simulation.  Future  applications  of  these  compressible 
subgrid-scale  models  to  the  large-eddy  simulation  of  more  complex  supersonic  flows  are  dis¬ 
cussed  briefly. 


Jacobs,  P.  A.,  R.  C.  Rogers,  E.  H.  Weidner,  and  R.  D.  Bittner:  Flow  establishment  in  a 
generic  scramjet  combustor.  ICASE  Report  No.  90-77,  October  16,  1990,  40  pages.  Sub¬ 
mitted  to  the  AIAA  Journal. 


The  establishment  of  a  quasi-steady  flow  in  a  generic  scramjet  combustor  is  studied  for 
the  case  of  a  time  varying  inflow  to  the  combustor.  Such  transient  flow  is  characteristic  of 
the  reflected-shock  tunnel  and  expansion-tube  test  facilities.  Several  numerical  simulations 
of  hypervelocity  flow  through  a  straight-duct  combusto*  /ith  either  a  side- wall-step  fuel 
injector  or  a  centrally-located  strut  injector  are  presented.  Comparisons  are  made  between 
impulsively  started  but  otherwise  constant  flow  conditions  (typical  of  the  expansion-tube 
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or  tailored  operation  of  the  reflected-shock  tunnel)  and  the  relaxing  flow  produced  by  the 
“undertailored”  operation  of  the  reflected-shock  tunnel.  Generally  the  inviscid  flow  features, 
such  as  the  shock  pattern  and  pressure  distribution,  were  unaffected  by  the  time  varying  inlet 
conditions  and  approached  steady  state  in  approximately  the  times  indicated  by  experimental 
correlations.  However,  viscous  features,  such  as  heat  transfer  and  skin  friction,  were  altered 
by  the  relaxing  inlet  flow  conditions. 

Papageorgiou,  Demetrius  T.,  and  Yiorgos  S.  Smyrlis:  The  route  to  chaos  for  the  Kuramoto- 
Sivashinsky  equation.  ICASE  Report  No.  90-78,  October  19,  1990,  55  pages.  Submitted  to 
Theoretical  and  Computational  Fluid  Dynamics. 


We  present  the  results  of  extensive  numerical  experiments  of  the  spatially  periodic  initial 
value  problem  for  the  Kuramoto-Sivashinsky  equation.  Our  concern  is  with  the  asymptotic 
nonlinear  dynamics  as  the  dissipation  parameter  decreases  and  spatio-temporal  chaos  sets 
in.  To  this  end  the  initial  condition  is  taken  to  be  the  same  for  all  numerical  experiments  (a 
single  sine  wave  is  used)  and  the  large  time  evolution  of  the  system  is  followed  numerically. 
Numerous  computations  were  performed  to  establish  the  existence  of  windows,  in  parameter 
space,  in  which  the  solution  has  the  following  characteristics  as  the  viscosity  is  decreased: 
A  steady  fully  modal  attractor  to  a  steady  bimodal  attractor  to  another  steady  fully  modal 
attractor  to  a  steady  trimodal  attractor  to  a  periodic  (in  time)  attractor,  to  another  steady 
fully  modal  attractor,  to  another  time-periodic  attractor,  to  a  steady  tetramodal  attractor,  to 
another  time-periodic  attractor  having  a  full  sequence  of  period-doublings  (in  the  parameter 
space)  to  chaos.  Numerous  solutions  are  presented  which  provide  conclusive  evidence  of 
the  periodic-doubling  cascades  which  precede  chaos  for  this  infinite-dimensional  dynamical 
system.  These  results  permit  a  computation  of  the  lengths  of  subwindows  which  in  turn 
provide  an  estimate  for  their  successive  ratios  as  the  cascade  develops.  A  calculation  based 
on  the  numerical  results  is  also  presented  to  show  that  the  period  doubling  sequences  found 
here  for  the  Kuramoto-Sivashinsky  equation,  are  in  complete  agreement  with  Feigenbaum’s 
universal  constant  of  4.669201609...  .  Some  preliminary  work  shows  several  other  windows 
following  the  first  chaotic  one  including  periodic,  chaotic  and  a  steady  octamodal  window; 
however  the  windows  shrink  significantly  in  size  to  enable  concrete  quantitative  conclusions 
to  be  made. 


Smith,  R.  C.,  K.  L.  Bowers  and  J.  Lund:  A  fully  Sinc-Galerkin  method  for  Euler- Bernoulli 
beam  models.  ICASE  Report  No.  90-79,  October  23,  1990,  44  pages.  Submitted  to  Numeri¬ 
cal  Methods  for  Partial  Differential  Equations. 


A  fully  Sinc-Galerkin  method  in  both  space  and  time  is  presented  for  fourth-order  time- 
dependent  partial  differential  equations  with  fixed  and  cantilever  boundary  conditions.  The 
sine  discretizations  for  the  second-order  temporal  problem  and  the  fourth-order  spatial  prob¬ 
lems  are  presented.  Alternate  formulations  for  variable  parameter  fourth-order  problems  are 
given  which  prove  to  be  especially  useful  when  applying  the  forward  techniques  of  this 
paper  to  parameter  recovery  problems.  The  discrete  system  which  corresponds  to  the  time- 
dependent  partial  differential  equations  of  interest  are  then  formulated.  Computational  issues 
are  discussed  and  a  robust  and  efficient  algorithm  for  solving  the  resulting  matrix  system  is 
outlined.  Numerical  results  which  highlight  the  method  are  given  for  problems  with  both 
analytic  and  singular  solutions  as  well  as  fixed  and  cantilever  boundary  conditions. 
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Piomelli,  Ugo,  and  Thomas  A.  Zang:  Large-eddy  simulation  of  transitional  channel  flow. 
ICASE  Report  No.  90-80,  November  5,  1990,  14  pages.  To  appear  in  Computer  Physics 
Communications. 


A  large-eddy  simulation  (LES)  of  transition  in  plane  channel  flow  has  been  carried  out. 
The  LES  results  have  been  compared  with  those  of  a  fine  direct  numerical  simulation  (DNS), 
and  with  those  of  a  coarse  DNS  that  uses  the  same  mesh  as  the  LES,  but  no  residual  stress 
model.  While  at  the  early  stages  of  transition  LES  and  coarse  DNS  give  the  same  results, 
the  presence  of  the  residual  stress  model  was  found  to  be  necessary  to  predict  accurately 
mean  velocity  and  Reynolds  stress  profiles  during  the  late  stages  of  transition  (after  the 
second  spike  stage).  The  evolution  of  single  Fourier  modes  is  also  predicted  more  accurately 
by  the  LES  than  by  the  DNS.  As  small  scales  are  generated,  the  dissipative  character  of 
the  residual  stress  starts  to  reproduce  correctly  the  energy  cascade;  as  transition  progresses, 
then,  and  the  flow  approaches  its  fully  developed  turbulent  state,  the  subgrid  scales  tend 
towards  equilibrium  and  the  model  becomes  more  accurate. 


Shaw,  Stephen  J.,  and  Peter  W.  Duck:  The  inviscid  stability  of  supersonic  flow  past  heated 
or  cooled  axisymmetric  bodies.  ICASE  Report  No.  90-81,  November  30,  1990,  53  pages. 
Submitted  to  Physics  of  Fluids. 


The  inviscid,  linear,  non-axisymmetric,  temporal  stability  of  the  boundary  layer  associ¬ 
ated  with  the  supersonic  flow  past  axisymmetric  bodies  (with  particular  emphasis  on  long 
thin,  straight  circular  cylinders),  subject  to  heated  or  cooled  wall  conditions  is  investigated. 
The  eigenvalue  problem  is  computed  in  some  detail  for  a  particular  Mach  number  of  3.8, 
revealing  that  the  effect  of  curvature  and  the  choice  of  wall  conditions  both  have  a  significant 
effect  on  the  stability  of  the  flow. 

Both  the  asymptotic,  large  azimuthal  wavenumber  solution  and  the  asymptotic,  far  down¬ 
stream  solution  are  obtained  for  the  stability  analysis  and  compared  with  numerical  results. 
Additionally,  asymptotic  analyses  valid  for  large  radii  of  curvature  with  cooled/heated  wall 
conditions,  are  presented.  We  find,  in  general,  important  differences  exist  between  the  wall 
temperature  conditions  imposed  in  this  paper  and  the  adiabatic  wall  conditions  considered 
previously. 


Jackson,  T.  L.,  and  C.  E.  Grosch:  Zero  wavenumber  modes  of  a  compressible  supersonic 
mixing  layer.  ICASE  Report  No.  90-82,  December  3,  1990,  12  pages.  To  be  submitted  to 
Physics  of  Fluids  A. 


In  this  paper  we  show  that  there  exists  a  family  of  supersonic  neutral  modes  for  a  com¬ 
pressible  mixing  layer  in  an  unbounded  domain.  These  modes  have  zero  wavenumber  and 
frequency  with  non-zero  phase  speed.  They  are  analogous  to  the  supersonic  neutral  modes 
of  the  compressible  vortex  sheet  found  by  Miles.  The  results  presented  here  give  a  more 
complete  picture  of  the  spectrum  of  the  disturbances  in  this  flow. 
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Liandrat,  J.,  and  Ph.  Tchamitchian:  Resolution  of  the  ID  regularized  Burgers  equation 
using  a  spatial  wavelet  approximation.  ICASE  Report  No.  90-83,  December  5,  1990,  36 
pages.  Submitted  to  Annales  Institut  Henri  Poincare. 


The  Burgers  equation  with  a  small  viscosity  term,  initial  and  periodic  boundary  condi¬ 
tions  is  resolved  numerically  using  a  spatial  approximation  constructed  from  an  orthonormal 
basis  of  wavelets. 

The  algorithm  is  directly  derived  from  th«  notions  of  multiresolution  analysis  and  tree 
algorithms.  Before  the  numerical  algorithm  is  described  these  notions  are  first  recalled. 
The  method  uses  extensively  the  localization  properties  of  the  wavelets  in  the  physical  and 
Fourier  spaces.  Moreover,  we  take  advantage  of  the  fact  that  the  involved  linear  operators 
have  constant  coefficients.  Finally,  the  algorithm  can  be  considered  as  a  time  marching 
version  of  the  tree  algorithm. 

The  most  important  point  is  that  an  adaptative  version  of  the  algorithm  exists:  it  allows 
one  to  reduce  in  a  significant  way  the  number  of  degrees  of  freedom  required  for  a  good 
computation  of  the  solution. 

Numerical  results  and  description  of  the  different  elements  of  the  algorithm  are  provided 
in  combination  with  different  mathematical  comments  on  the  method  and  some  comparison 
with  more  classical  numerical  algorithms. 


Denier,  J.  P.,  and  Philip  Hall:  The  effect  of  wall  compliance  on  the  Gortler  vortex  instability. 
ICASE  Report  No.  90-84,  December  6,  1990,  15  pages.  Submitted  to  Physics  of  Fluids. 


The  stability  of  the  flow  of  a  viscous  incompressible  fluid  over  a  curved  compliant  wall  to 
longitudinal  Gortler  vortices  is  investigated.  The  compliant  wail  is  modeled  by  a  particularly 
simple  equation  relating  the  induced  wall  displacement  to  the  pressure  in  the  overlying  fluid. 
Attention  is  restricted  to  the  large  Gortler  number  regime;  this  regime  being  appropriate  to 
the  most  unstable  Gortler  mode.  The  effect  of  wall  compliance  on  this  most  unstable  mode 
is  investigated. 


Fu,  Yibin  B.,  Philip  Hall,  and  Nicholas  D.  Blackaby:  On  the  Gortler  instability  in  hypersonic 
flows:  Sutherland  law  fluids  and  real  gas  effects.  ICASE  Report  No.  90-85,  December  11, 
1990,  70  pages.  Submitted  to  Proc.  Roy.  Soc. 


The  Gortler  vortex  instability  mechanism  in  a  hypersonic  boundary  layer  on  a  curved  wall 
is  investigated  in  this  paper.  Our  aim  is  to  clarify  the  precise  roles  of  the  effects  of  boundary 
layer  growth,  wall  cooling  and  gas  dissociation  in  the  determination  of  stability  properties. 
We  first  assume  that  the  fluid  is  an  ideal  gas  with  viscosity  given  by  Sutherland’s  law.  It  is 
shown  that  when  the  free  stream  Mach  number  M  is  large,  the  boundary  layer  divides  into 
two  sublayers:  a  wall  layer  of  0(M3/2)  thickness  over  which  the  basic  state  temperature  is 
0(M2)  and  a  temperature  adjustment  layer  of  0(1)  thickness  over  which  the  basic  tempera¬ 
ture  decreases  monotonically  to  its  free  stream  value.  Gortler  vortices  which  have  wavelength 
comparable  with  the  boundary  layer  thickness  (i.e.,  have  local  wavenumber  of  order  M-3/2) 
are  referred  to  as  wall  modes.  We  show  that  their  downstream  evolution  is  governed  by  a  set 
of  parabolic  partial  differential  equations  and  that  they  have  the  usual  features  of  Gortler 
vortices  in  incompressible  boundary  layers.  As  the  local  wavenumber  increases,  the  neutral 
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Gortler  number  decreases  and  the  centre  ofvortex  activity  moves  towards  the  temperature 
adjustment  layer.  Gortler  vortices  with  wavenumber  of  order  one  or  larger  must  necessarily 
be  trapped  in  the  temperature  adjustment  layer  and  it  is  this  mode  which  is  the  most  dan¬ 
gerous.  For  this  mode,  we  find  that  the  leading  order  term  in  the  Gortler  number  expansion 
is  independent  of  the  wavenumber  and  is  due  to  the  curvature  of  the  basic  state.  This  term 
is  also  the  asymptotic  limit  of  the  neutral  Gortler  numbers  of  the  wall  mode.  To  determine 
the  higher  ordercorrection  terms  in  the  Gortler  number  expansion,  we  have  to  distinguish 
between  two  wall  curvature  cas^;.  When  the  wall  curvature  is  proportional  to  (2a:)-3/2  where 
x  is  the  streamwise  variable,  the  Mach  number  M  can  be  scaled  out  of  the  problem  and  we 
show  that  in  the  0(1)  wavenumber  regime,  Gortler  vortices  are  again  governed  by  a  set  of 
parabolicpartial  differential  equations  and  therefore  the  higher  order  correction  terms  in  the 
Gortler  number  expansion  are  not  uniquely  determined  and  are  strongly  dependent  on  non- 
parallel  effects.  In  the  large  wavenumber  limit,  however,  nonparallel  effects  become  of  second 
order;  Gortler  vortices  evolve  downstream  in  a  quasiparallel  manner  and  the  Gortler  number 
expansion  has  its  first  three  terms  independent  of  nonparallel  effects.  In  the  more  general 
case  when  the  wall  curvature  is  not  proportional  to  (2x)~3//2,  the  effect  of  the  curvature 
of  the  basic  state  persists  in  the  downstream  development  of  Gortler  vortices;  non-parallel 
effects  are  important  over  a  larger  range  of  wavenumbers  and  they  become  a  second  order 
only  when  the  wavenumber  is  of  order  higher  than  0(M 1/f4).  In  the  latter  case  the  Gortler 
number  expansion  has  the  first  two  terms  independent  of  nonparallel  effects;  the  first  term 
being  due  to  the  curvature  of  the  basic  state  and  the  second  term  due  to  viscous  effects. 
The  second  term  becomes  comparable  with  the  first  term  when  the  wavenumber  reaches 
the  order  M3/8,  in  which  case  another  correction  term  can  also  be  found  independently  of 
nonparallel  effects.  Next  we  investigate  real  gas  effects  by  assuming  that  the  fluid  is  an  ideal 
dissociating  gas.  We  find  that  both  gas  dissociation  and  wall  cooling  are  destabilizing  for 
the  mode  trapped  in  the  temperature  adjustment  layer,  but  for  the  wall  mode  trapped  near 
the  wall  the  effect  of  gas  dissociation  can  be  either  destabilizing  or  stabilizing. 


Hall,  Philip,  and  P.  Balakumar:  On  a  class  of  unsteady  three-dimensional  Navier  Stokes 
solutions  relevant  to  rotating  disc  flows:  Threshold  amplitudes  and  finite  time  singularities. 
ICASE  Report  No.  90-86  December  10,  1990,  45  pages.  Submitted  to  JFM. 


A  class  of  exact  steady  and  unsteady  solutions  of  the  Navier  Stokes  equations  in  cylindri¬ 
cal  polar  coordinates  is  given.  The  flows  correspond  to  the  motion  induced  by  an  infinite  disc 
rotating  with  constant  angular  velocity  about  the  z-axis  in  a  fluid  occupying  a  semi-infinite 
region  which,  at  large  distances  from  the  disc,  has  velocity  field  proportional  to  (x,  —y,  0) 
with  respect  to  a  Cartesian  coordinate  system.  It  is  shown  that  when  the  rate  of  rotation 
is  large  Karman’s  exact  solution  for  a  disc  rotating  in  an  otherwise  motionless  fluid  is  re¬ 
covered.  In  the  limit  of  zero  rotation  rate  a  particular  form  of  Howarth’s  exact  solution 
for  three-dimensional  stagnation  point  flow  is  obtained.  The  unsteady  form  of  the  partial 
differential  system  describing  this  class  of  flow  may  be  generalized  to  time-periodic  equilib¬ 
rium  flows.  In  addition  the  unsteady  equations  are  shown  to  describe  a  strongly  nonlinear 
instability  of  Karman’s  rotating  disc  flow.  It  is  shown  that  sufficiently  large  perturbations 
lead  to  a  finite  time  breakdown  of  that  flow  whilst  smaller  disturbances  decay  to  zero.  If  the 
stagnation  point  flow  at  infinity  is  sufficiently  strong  the  steady  basic  states  become  linearly 
unstable.  In  fact  there  is  then  a  continuous  spectrum  of  unstable  eigenvalues  of  the  stability 
equations  but,  if  the  initial  value  problem  is  considered,  it  is  found  that,  at  large  values  of 
time,  the  continuous  spectrum  leads  to  a  velocity  field  growing  exponentially  in  time  with 
an  amplitude  decaying  algebraically  in  time. 
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Garbey,  Marc,  and  Jeffrey  S.  Scroggs:  Asymptotic-Induced  numerical  methods  for  conserva¬ 
tion  laws.  ICASE  Report  No.  90-87,  December  12,  1990,  24  pages.  Proceedings  of  Asymp¬ 
totic  Analysis  and  Numerical  Solution  of  Partial  Differential  Equations. 


Asymptotic-induced  methods  are  presented  for  the  numerical  solution  of  hyperbolic  con¬ 
servation  laws  with  or  without  viscosity.  The  methods  consist  of  multiple  stages.  The  first 
stage  is  to  obtain  a  first  approximation  by  using  a  first-order  method,  such  as  the  Godunov 
scheme.  Subsequent  stages  of  the  method  involve  solving  internal-layer  problems  identified 
by  using  techniques  derived  via  asymptotics.  Finally,  a  residual  correction  increases  the 
accuracy  of  the  scheme.  The  method  is  derived  and  justified  with  singular  perturbation 
techniques. 


Speziale,  Charles  G.,  Rishi  Raj,  and  Thomas  B.  Gatski:  Modeling  the  dissipation  rate  in 
rotating  turbulent  flows.  ICASE  Report  No.  90-88,  December  12,  1990,  28  pages.  To  appear 
in  Recent  Developments  in  Turbulence,  (Springer- Verlag). 


A  variety  of  modifications  to  the  modeled  dissipation  rate  transport  equation  that  have 
been  proposed  during  the  past  two  decades  to  account  for  rotational  strains  are  examined. 
The  models  are  subjected  to  two  crucial  test  cases:  the  decay  to  isotropic  turbulence  in  a 
rotating  frame  and  homogeneous  shear  flow  in  a  rotating  frame.  It  is  demonstrated  that 
these  modifications  do  not  yield  substantially  improved  predictions  for  these  two  test  cases 
and  in  many  instances  give  rise  to  unphysical  behavior.  An  alternative  proposal,  based  on 
the  use  of  the  tensor  dissipation  rate,  is  made  for  the  development  of  improved  models. 


Balakumar,  P.,  P.  Hall,  and  M.  R.  Malik:  On  the  receptivity  and  non-parallel  stability  of 
travelling  disturbances  in  rotating  disk  flow.  ICASE  Report  No.  90-89,  December  12,  1990, 
31  pages.  Submitted  to  Theoretical  and  Computational  Fluid  Dynamics. 


The  generation  and  evolution  of  small  amplitude  long  wavelength  travelling  disturbances 
in  rotating  disk  flow  is  the  subject  of  this  paper.  The  steady  rotational  speed  of  the  disk 
is  perturbed  so  as  to  introduce  high  frequency  oscillations  in  the  flow  field.  Secondly,  we 
introduce  surface  imperfections  on  the  disk  such  as  roughness  elements.  The  interaction 
of  these  two  disturbances  will  generate  the  instabilty  waves  whose  evolution  is  governed 
by  parabolic  partial  differential  equations  which  are  solved  numerically.  It  is  found  that, 
for  the  class  of  disturbances  considered  here  (wavelength  on  the  order  of  Reynolds  number), 
eigensolutions  exist  which  decay  or  grow  algebraically  in  the  radial  direction.  However,  these 
solutions  grow  only  for  frequencies  larger  than  4.58  times  the  steady  rotational  speed  of  the 
disk.  The  computed  receptivity  coefficient  shows  that  there  is  an  optimum  size  of  roughness 
for  which  these  modes  are  excited  the  most.  The  width  of  these  roughness  elements  in  the 
radial  direction  is  about  .lrjj  where  rj  is  the  radial  location  of  the  roughness.  It  is  also 
found  that  the  receptivity  coefficient  is  larger  for  a  negative  spanwise  wavenumber  than  for 
a  positive  one.  Typical  wave  angles  found  for  these  disturbances  are  about  —26°. 
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Narasimha,  R.:  Modelling  the  transitional  boundary  layer.  ICASE  Report  No.  90-90)  De¬ 
cember  17,  1990,  22  pages.  Submitted  to  Applied  Mechanics  Review. 


Recent  developments  in  the  modelling  of  the  transition  zone  in  the  boundary  layer  are 
reviewed  (the  zone  being  defined  as  extending  from  the  station  where  intermittency  begins  to 
depart  from  zero  to  that  where  it  is  nearly  unity).  The  value  of  using  a  new  non-dimensional 
spot  formation  rate  parameter,  and  the  importance  of  allowing  for  so-called  subtransitions 
within  the  transition  zone,  are  both  stressed.  Models  do  reasonably  wel)  in  constant  pressure 
2-dimensional  flows,  but  in  the  presence  of  strong  pressure  gradients  further  improvements 
are  needed.  The  linear  combination  approach  works  surprisingly  well  in  most  cases,  but 
would  not  be  so  successful  in  situations  where  a  purely  laminar  boundary  layer  would  separate 
but  a  transitional  one  would  not.  Intermittency-weighted  eddy  viscosity  methods  ao  not 
predict  peak  surface  parameters  well  without  the  introduction  of  an  overshooting  transition 
function  whose  connection  with  the  spot  theory  of  transition  is  obscure. 

Suggestions  are  made  for  further  work  that  now  appears  necessary  for  developing  im¬ 
proved  models  of  the  transition  zone. 


Jacobs,  Peter  A.:  Transient,  hypervelocity  flow  in  an  axisymmetric  nozzle.  ICASE  Report 
No.  91-1.  January  2,  1991,  52  pages.  To  be  submitted  as  an  AIAA  paper. 


This  study  examines  the  performance  of  an  axisymmetric  nozzle  which  was  designed 
to  produce  uniform,  parallel  flow  with  a  nominal  Mach  number  of  8.  A  free-piston-driven 
shock  tube  was  used  to  supply  the  nozzle  with  high-temperature,  high-pressure  test  gas. 
Performance  was  assessed  by  measuring  Pitot  pressures  across  the  exit  plane  of  the  nozzle 
and,  over  the  range  of  operating  conditions  examined,  the  nozzle  produced  satisfactorily  test 
flows.  However,  there  were  flow  disturbances  that  persisted  for  significant  times  after  flow 
initiation. 

The  detailed  starting  process  of  the  nozzle  was  also  investigated  by  performing  numeri¬ 
cal  simulations  at  several  nominal  test  conditions.  The  classical  description  of  the  starting 
process,  based  on  a  quasi-one-dimensional  model,  provided  a  reasonable  approximation  and 
was  used  to  demonstrate  that  the  starting  process  could  consume  a  significant  fraction  of 
the  otherwise  usable  test  gas.  This  was  especially  important  at  high  operating  enthalpies 
where  nozzle  supply  conditions  were  maintained  for  shorter  times.  Multidimensional  simu¬ 
lations  illustrated  a  mechanism  by  which  the  starting  process  in  the  actual  nozzle  could  take 
longer  than  that  predicted  by  the  quasi-one-dimensional  analysis.  However,  the  cause  of  the 
persistent  disturbances  observed  in  the  experimental  calibration  was  not  identified. 


Ta’asan,  Shlomo:  “One  shot  ”  methods  for  optimal  control  of  distributed  parameters  systems  I: 
Finite  dimensional  control.  ICASE  Report  No.  91-2,  January  8,  1991,  22  pages.  Submitted 
to  SIAM  Journal  on  Optimization. 


This  paper  discusses  the  efficient  numerical  treatment  of  optimal  control  problems  gov¬ 
erned  by  elliptic  PDE’s  and  systems  of  elliptic  PDE’s,  where  the  control  is  finite  dimensional. 
Distributed  control  as  well  as  boundary  control  cases  are  discussed.  The  main  characteristic 
of  the  new  methods  is  that  they  are  designed  to  solve  the  full  optimization  problem  directly, 
rather  than  accelerating  a  descent  method  by  an  efficient  multigrid  solver  for  the  equations 
involved.  The  methods  use  the  adjoint  state  in  order  to  achieve  efficient  smoother  and  a 
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robust  coarsening  strategy.  The  main  idea  is  the  treatment  of  the  control  variables  on  ap¬ 
propriate  scales,  i.e.,  control  variables  that  correspond  to  smooth  functions  are  solved  for  on 
coarse  grids  depending  on  the  smoothness  of  these  functions.  Solution  of  the  control  prob¬ 
lems  is  achieved  with  the  cost  of  solving  the  constraint  equations  about  two  to  three  times 
(by  a  multigrid  solver).  Numerical  examples  demonstrate  the  effectiveness  of  the  method 
proposed  in  distributed  control  case,  pointwise  control  and  boundary  control  problems. 


Crockett,  Thomas  W.,  and  Tobias  Orloff:  A  parallel  rendering  algorithm  for  MIMD  archi¬ 
tectures.  ICASE  Report  No.  91-3. 


Applications  such  as  animation  and  scientific  visualization  demand  high  performance  ren¬ 
dering  of  complex  three  dimensional  scenes.  To  deliver  the  necessary  rendering  rates,  highly 
parallel  hardware  architectures  are  required.  The  challenge  is  then  to  design  algorithms  and 
software  which  effectively  use  the  hardware  parallelism.  This  paper  describes  a  rendering 
algorithm  targeted  to  distributed  memory  MIMD  architectures.  For  maximum  performance, 
the  algorithm  exploits  both  object-level  and  pixel-level  parallelism.  The  behavior  of  the  algo¬ 
rithm  is  examined  both  analytically  and  experimentally.  Its  performance  for  large  numbers 
of  processors  is  found  to  be  limited  primarily  by  communication  overheads.  An  experimental 
implementation  for  the  Intel  iPSC/860  shows  increasing  performance  from  1  to  128  proces¬ 
sors  across  a  wide  range  of  scene  complexities.  It  is  shown  that  minimal  modifications  to 
the  algorithm  will  adapt  it  for  use  on  shared  memory  architectures  as  well. 


Bokhari,  Shahid  H.:  Complete  exchange  on  the  iPSC-860.  ICASE  Report  No.  91-4,  January 
11,  1991,  32  pages.  Submitted  to  the  1991  International  Conference  on  Parallel  Processing. 


The  implementation  of  complete  exchange  on  the  circuit  switched  Intel  iPSC-860  hyper¬ 
cube  is  described.  This  pattern,  also  known  as  all-to-all  personalized  communication,  is  the 
densest  requirement  that  can  be  imposed  on  a  network.  On  the  iPSC-860,  care  needs  to  be 
taken  to  avoid  edge  contention,  which  can  have  a  disastrous  impact  on  communication  time. 
There  are  basically  two  classes  of  algorithms  that  achieve  contention  free  complete  exchange. 
The  first  contains  the  classical  standard  exchange  algorithm  that  is  generally  useful  for  small 
message  sizes.  The  second  includes  a  number  of  optimal  or  near-optimal  algorithms  that  are 
best  for  large  messages. 

Measurements  of  communication  overhead  on  the  iPSC-860  are  given  and  a  notation 
for  analyzing  communication  link  usage  is  developed.  It  is  shown  that  for  the  two  classes 
of  algorithms,  there  is  substantial  variation  in  performance  with  synchronization  technique 
and  choice  of  message  protocol.  Timings  of  six  implementations  are  given;  each  of  these  is 
useful  over  a  particular  range  of  message  size  and  cube  dimension. 

Since  the  complete  exchange  is  a  superset  of  all  communication  patterns,  these  timings 
represent  upper  bounds  on  the  time  required  by  an  arbitrary  communication  requirement. 
These  results  indicate  that  the  programmer  needs  to  evaluate  several  possibilities  before 
finalizing  an  implementation — a  careful  choice  can  lead  to  very  significant  savings  in  time. 
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Bokhari,  Shahid  H.:  Multiphase  complete  exchange  on  a  circuit  switched  hypercube.  ICASE 
Report  No.  91-5,  January  8,  1991,  24  pages.  Submitted  to  1991  International  Conference  on 
Parallel  Processing. 

On  a  distributed  memory  parallel  computer,  the  complete  exchange  (all-to-all  personal¬ 
ized)  communication  pattern  requires  each  of  n  processors  to  send  a  different  block  of  data 
to  each  of  the  remaining  n  —  1  processors.  This  pattern  is  at  the  heart  of  many  important 
algorithms,  most  notably  the  matrix  transpose. 

For  a  circuit  switched  hypercube  of  dimension  d  (n  =  2d),  two  algorithms  for  achieving 
complete  exchange  are  known.  These  are  (1)  the  Standard  Exchange  approach  that  employs 
d  transmissions  of  size  2d~l  blocks  each  and  is  useful  for  small  block  sizes,  and  (2)  the  Optimal 
Circuit  Switched  algorithm  that  employs  2^  —  1  transmissions  of  1  block  each  and  is  best  for 
large  block  sizes. 

A  unified  multiphase  algorithm  is  described  that  includes  these  two  algorithms  as  special 
cases.  The  complete  exchange  on  a  hypercube  of  dimension  d  and  block  size  m  is  achieved 
by  carrying  out  k  partial  exchanges  on  subcubes  of  dimension  d{,  Ej !=1d;  =  d  and  effective 
block  size  m<  =  m2d~di.  When  k  =  d  and  all  d,  =  1,  this  corresponds  to  algorithm  (1)  above. 
For  the  case  of  k  =  1  and  d\  =  d,  this  becomes  the  circuit  switched  algorithm  (2).  Changing 
the  subcube  dimensions  d,  varies  the  effective  block  size  and  permits  a  compromise  between 
the  data  permutation  and  block  transmission  overhead  of  (1)  and  the  startup  overhead  of 
(2). 

For  a  hypercube  of  dimension  d,  the  number  of  possible  combinations  of  subcubes  is  p(d), 
the  number  of  partitions  of  the  integer  d.  This  is  an  exponential  but  very  slowly  growing 
function  (e.g.  p( 7)  =  15,p(10)  =  42)  and  it  is  feasible  to  enumerate  over  these  partitions  to 
discover  the  best  combination  for  a  given  message  size. 

This  approach  has  been  analyzed  for,  and  implemented  on,  the  Intel  iPSC-860  circuit 
switched  hypercube.  Measurements  show  good  agreement  with  predictions  and  demonstrate 
that  the  multiphase  approach  can  substantially  improve  performance  for  block  sizes  in  the  0- 
160  byte  range.  This  range,  which  corresponds  to  0-40  floating  point  numbers  per  processor, 
is  commonly  encountered  in  practical  numeric  applications.  The  multiphase  technique  is 
applicable  to  all  circuit-switched  hypercubes  that  use  the  common  1  e-cube’  routing  strategy. 


Duck,  Peter  W.,  and  Mehdi  R.  Khorrami:  On  the  effects  of  viscosity  on  the  stability  of  a 
trailing-line  vortex.  ICASE  Report  No.  91-6,  January  10,  1991,  32  pages.  Submitted  to  J. 
Fluid  Mechanics. 


The  linear  stability  of  the  Batchelor  (1964)  voitex  is  investigated.  Particular  emphasis 
is  placed  on  modes  found  recently  in  a  numerical  study  by  Khorrami  (1991).  These  modes 
have  a  number  of  features  very  distinct  from  those  found  previously  for  this  vortex,  in¬ 
cluding  (i)  exhibiting  small  growth  rates  at  large  Reynolds  numbers  and  (ii)  susceptibility 
to  desiabilisation  by  viscosity.  In  this  paper  these  modes  are  described  using  asymptotic 
techniques,  producing  results  which  compare  very  favourably  with  fully  numerical  results  at 
large  Reynolds  numbers. 
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Jordan,  Harry  F.:  Shared  versus  distributed  memory  multiprocessors.  ICASE  Report  No. 
91-7.  January  10,  1991,  19  pages.  To  appear  in  the  Proceedings  of  European  Centre  for 
Medium  Range  Weather  Forecasts  workshop  on  Use  of  Parallel  Processors  in  Meteorology, 
Nov.  26-30,  1990. 


The  question  of  whether  multiprocessors  should  have  shared  or  distributed  memory  has 
attracted  a  great  deal  of  attention.  Some  researchers  argue  strongly  for  building  distributed 
memory  machines,  while  others  argue  just  as  strongly  for  programming  shared  memory  mul¬ 
tiprocessors.  A  great  deal  of  research  is  underway  on  both  types  of  parallel  systems.  This 
paper  puts  special  emphasis  on  systems  with  a  very  large  number  of  processsors  for  compu¬ 
tation  intensive  tasks  and  considers  research  and  implementation  trends.  It  appears  that  the 
two  types  of  system  will  likely  converge  to  a  common  form  for  large  scale  multiprocessors. 


Harden,  Ami:  Recent  developments  in  shock-capturing  schemes.  ICASE  Report  No.  91-8, 
January  22, 1991,  13  pages.  To  appear  in  Proc.  of  the  International  Congress  of  Mathemati¬ 
cians  1990,  Kyoto,  Japan. 


In  this  paper  we  review  the  development  of  the  shock-capturing  methodology,  paying 
special  attention  to  the  increasing  nonlinearity  in  its  design  and  its  relation  to  interpolation. 
It  is  well-known  that  high-order  approximations  to  a  discontinuous  function  generate  spurious 
oscillations  near  the  discontinuity  (Gibbs  phenomenon).  Unlike  standard  finite-difference 
methods  which  use  a  fixed  stencil,  modern  shock-capturing  schemes  use  an  adaptive  stencil 
which  is  selected  according  to  the  local  smoothness  of  the  solution.  Near  discontinuities 
this  technique  automatically  switches  to  one-sided  approximations,  thus  avoiding  the  use  of 
discontinuous  data  which  brings  about  spurious  oscillations. 


Speziale,  Charles  G.,  and  Sutanu  Sarkar:  Second-order  closure  models  for  supersonic  tur¬ 
bulent  flows.  ICASE  Report  No.  91-9,  January  23,  1991,  24  pages.  AIAA  29th  Aerospace 
Sciences  Meeting. 


Recent  work  by  the  authors  and  their  colleagues  on  the  development  of  a  second-order 
closure  model  for  high-speed  compressible  flows  is  reviewed.  This  turbulence  closure  is  based 
on  the  solution  of  modeled  transport  equations  for  the  Favre-averaged  Reynolds  stress  tensor 
and  the  solenoidal  part  of  the  turbulent  dissipation  rate.  A  new  model  for  the  compressible 
dissipation  is  used  along  with  traditional  gradient  transport  models  for  the  Reynolds  heat  flux 
and  mass  flux  terms.  Consistent  with  simple  asymptotic  analyses,  the  deviatoric  part  of  the 
remaining  higher-order  correlations  in  the  Reynolds  stress  transport  equation  are  modeled 
by  a  variable  density  extension  of  the  newest  incompressible  models.  The  resulting  second- 
order  closure  model  is  tested  in  a  variety  of  compressible  turbulent  flows  which  include  the 
decay  of  isotropic  turbulence,  homogeneous  shear  flow,  the  supersonic  mixing  layer,  and  the 
supersonic  flat-plate  turbulent  boundary  layer.  Comparisons  between  the  model  predictions 
and  the  results  of  physical  and  numerical  experiments  are  quite  encouraging. 
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Chen,  Yen-Ming,  Yuh-Roung  Ou,  and  Arne  J.  Pearlstein:  Development  of  the  wake  behind 
a  circular  cylinder  impulsively  started  into  rotatory  and  rectilinear  motion:  Intermediate 
rotation  rates.  ICASE  Report  No.  91-10,  January  24,  1991,  61  pages.  Submitted  to  Journal 
of  Fluid  Mechanics. 


The  temporal  development  of  two-dimensional  viscous  incompressible  flow  generated  by  a 
circular  cylinder  started  impulsively  into  steady  rotatory  and  rectilinear  motion  is  studied  by 
integration  of  a  velocity/vorticity  formulation  of  the  governing  equations,  using  an  explicit 
finite-difference/pseudo-spectral  technique  and  a  new  implementation  of  the  Biot-Savart  law. 
Results  axe  presented  for  a  Reynolds  number  of  200  (based  on  the  cylinder  diameter  2a  and 
the  magnitude  U  of  the  rectilinear  velocity)  for  several  values  of  the  angular/rectilinear 
speed  ratio  a  =  Clu/U  (where  ft  is  the  angular  speed)  up  to  3.25.  For  values  of  a  >  1, 
out  extension  of  the  computations  to  larger  dimensionless  times  than  those  possible  in  the 
experimental  work  of  Badr  &  Dennis  (1985)  allows  for  a  more  complete  discussion  of  the 
long-term  development  of  the  wake.  We  also  discuss  several  aspects  of  the  kinematics  and 
dynamics  of  the  flow  not  considered  earlier.  For  higher  values  of  a,  our  results  indicate  that 
for  Re  =  200,  vortex  shedding  does  indeed  occur  for  a  =  3.25  (and  possibly  for  higher  values 
of  a  also),  in  contrast  to  the  conclusion  of  Coutanceau  &  M’enard  (1985).  The  shedding 
process  is,  however,  very  different  from  that  which  gives  rise  to  the  usual  K’arm’an  vortex 
street  for  a  =  0.  In  particular,  consecutive  vortices  shed  by  the  body  can  be  shed  from  the 
same  side,  and  be  of  the  same  sense,  in  contrast  to  the  nonrotating  case,  in  which  mirror- 
image  vortices  of  opposite  sense  are  shed  alternately  on  opposite  sides  of  the  body.  The 
implications  of  the  results  are  discussed  in  relation  to  the  possibility  of  suppressing  vortex 
shedding  by  open-  or  closed-loop  control  of  the  rotation  rate. 


Mavriplis,  D.  J.,  and  L.  Martinelli:  Multigrid  solution  of  compressible  turbulent  flow  on 
unstructured  meshes  using  a  two-equation  model  ICASE  Report  No.  91-11,  January  31, 
1991,  41  pages.  Submitted  to  AIAA  Journal. 


The  steady-state  solution  of  the  system  of  equations  consisting  of  the  full  Navier-Stokes 
equations  and  two  turbulence  equations  has  been  obtained  using  a  multigrid  strategy  on 
unstructured  meshes.  The  flow  equations  and  turbulence  equations  are  solved  in  a  loosley 
coupled  manner.  The  flow  equations  are  advanced  in  time  using  a  multi-stage  Runge-Kutta 
time  stepping  scheme  with  a  stability  bound  local  time-step,  while  the  turbulence  equa¬ 
tions  are  advanced  in  a  point-implicit  scheme  with  a  time-step  which  guarantees  stability 
and  positivity.  Low  Reynolds  number  modifications  to  the  original  two-equation  model  are 
incorporated  in  a  manner  which  results  in  well  behaved  equations  of  arbitrarily  small  wall 
distances.  A  variety  of  aerodynamic  flows  are  solved  for,  initializing  all  quantities  with  uni¬ 
form  freestream  values.  Rapid  and  uniform  convergence  rates  for  the  flow  and  turbulence 
equations  are  observed. 


Saltz,  Joel,  Serge  Petiton,  and  Harry  Berryman:  Performance  effects  of  irregular  commu¬ 
nications  patterns  on  massively  parallel  multiprocessors.  ICASE  Report  No.  91-12,  January 
31,  1991,  22  pages.  Journal  of  Parallel  and  Distributed  Computing. 


We  conduct  a  detailed  study  of  the  performance  effects  of  irregular  communications 
patterns  on  the  CM-2.  We  characterize  the  communications  capabilities  of  the  CM-2  under 
a  variety  of  controlled  conditions. 


39 


In  the  process  of  carrying  out  our  performance  evaluation,  we  develop  and  make  extensive 
use  of  a  parameterized  synthetic  mesh.  In  addition  we  carry  out  timings  with  unstructured 
meshes  generated  for  aerodynamic  codes  and  a  set  of  sparse  matrices  with  banded  patterns 
of  non-zeros.  This  benchmarking  suite  stresses  the  communications  capabilities  of  the  CM-2 
in  a  range  of  different  ways.  Our  benchmark  results  demonstrate  that  it  is  possible  to  make 
effective  use  of  much  of  the  massive  concurrency  available  in  the  communications  network. 


Wu,  Janet,  Joel  Saltz,  Harry  Berryman,  and  Seema  Hiranandani:  Distributed  memory  com¬ 
piler  design  for  sparse  problems.  ICASE  Report  No.  91-13,  January  31,  1991,  40  pages. 
Submitted  to  IEEE  Trans.  Software  Engng. 


In  this  paper  we  describe  and  demonstrate  a  compiler  and  runtime  support  mechanism. 
The  methods  presented  here  are  capable  of  solving  a  wide  range  of  sparse  and  unstructured 
problems  in  scientific  computing.  The  compiler  takes  as  input  a  Fortran  77  program  enhanced 
with  specifications  for  distributing  data,  and  the  compiler  outputs  a  message  passing  program 
that  runs  on  a  distributed  memory  computer.  The  runtime  support  for  this  compiler  is  a 
library  of  primitives  designed  to  efficiently  support  irregular  patterns  of  distributed  array 
accesses  and  irregular  distributed  array  partitions.  We  present  a  variety  of  Intel  iPSC/860 
performance  results  obtained  through  the  use  of  this  compiler. 


Smith,  R.C.,  K.L.  Bowers,  and  C.R.  Vogel:  Numerical  recovery  of  material  parameters  in 
Euler-Bernoulli  beam  models.  ICASE  Report  No.  91-14,  February  5,  1991,  38  pages.  Sub¬ 
mitted  to  Journal  of  Mathematical  Systems  and  Control. 


A  fully  Sine- Galer kin  method  for  recovering  the  spatially  varying  stiffness  parameter 
in  fourth-order  time-dependent  problems  with  fixed  and  cantilever  boundary  conditions  is 
presented.  The  forward  problems  are  discretized  with  a  sine  basis  in  both  the  spatial  and 
temporal  domains.  This  yields  an  approximate  solution  which  converges  exponentially  and 
is  valid  on  the  infinite  time  interval.  When  the  forward  methods  are  applied  to  parameter 
recovery  problems,  the  resulting  inverse  problems  are  ill-posed.  Tikhonov  regularization  is 
applied  and  the  resulting  minimization  problems  are  solved  via  a  quasi-Newton/trust  region 
algorithm.  The  L-curve  method  is  used  to  determine  an  appropriate  value  of  the  regulariza¬ 
tion  parameter.  Numerical  results  which  highlight  the  method  are  given  for  problems  with 
both  fixed  and  cantilever  boundary  conditions. 


Smith,  R.C.,  and  K.L.  Bowers:  A  fully  Galerkin  method  for  the  recovery  of  stiffness  and 
damping  parameters  in  Euler-Bernoulli  beam  models.  ICASE  Report  No.  91-15,  February 
5,  1991,  22  pages.  Submitted  to  Proc.  of  the  Second  Conference  on  Computation  and 
Control. 


A  fully  Sinc-Galerkin  method  for  recovering  the  spatially  varying  stiffness  and  damp¬ 
ing  parameters  in  Euler-Bernoulli  beam  models  is  presented.  The  forward  problems  are 
discretized  with  a  sine  basis  in  both  the  spatial  and  temporal  domains  thus  yielding  an  ap¬ 
proximate  solution  which  converges  exponentially  and  is  valid  on  the  infinite  time  interval. 
Hence  the  method  avoids  the  time-stepping  which  is  characteristic  of  many  of  the  forward 
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schemes  which  are  employed  in  parameter  recovery  algorithms.  Tikhonov  regularization  is 
used  to  stabilize  the  resulting  inverse  problem,  and  the  L-curve  method  for  determining 
an  appropriate  value  of  the  regularization  parameter  is  briefly  discussed.  Numerical  exam¬ 
ples  are  given  which  demonstrate  the  applicability  of  the  method  for  both  individual  and 
simultaneous  recovery  of  the  material  parameters. 


Thangam,  S.,  R.  Abid,  and  C.G.  Speziale:  Application  of  a  new  K  —  r  model  to  near  wall 
turbulent  flows.  ICASE  Report  No.  91-16,  February  5,  1991,  11  pages.  To  be  submitted  to 
AIAA  Journal. 


A  recently  developed  K  —  t  model  for  near  wall  turbulent  flows  is  applied  to  two  severe 
test  cases.  The  turbulent  flows  considered  include  the  incompressible  flat  plate  bound¬ 
ary  layer  with  adverse  pressure  gradients  and  incompressible  flow  past  a  backward  facing 
step.  Calculations  are  performed  for  this  two-equation  model  using  an  anisotropic  as  well  as 
isotropic  eddy-viscosity.  The  model  predictions  are  shown  to  compare  quite  favourably  with 
experimental  data. 


Balachandar,  S.,  and  L.  Sirovich:  Probability  distribution  functions  in  turbulent  convection. 
ICASE  Report  No.  91-17,  February  11,  1991,  40  pages.  Submitted  to  Physics  of  Fluids. 


Results  of  an  extensive  investigation  of  probability  distribution  functions  (pdfs)  for 
Rayleigh-Benard  convection,  in  the  hard  turbulence  regime,  is  presented.  It  is  seen  that 
the  pdfs  exhibit  a  high  degree  of  internal  universality.  In  certain  cases  this  universality  is 
established  within  two  Kolmogorov  scales  of  a  boundary.  A  discussion  of  the  factors  leading 
to  universality  is  presented. 


Jordan,  Harry  F.:  Digital  optical  computers  at  the  optoelectric  computing  systems  center. 
ICASE  Report  No.  91-18,  February  11,  1991,  20  pages.  Submitted  to  the  International 
Congress  on  Optical  Science  and  Engineering,  ECO’4,  The  Hague,  The  Netherlands,  March 
11-15,  1991. 


The  Digital  Optical  Computing  Program  within  the  National  Science  Foundation  En¬ 
gineering  Research  Center  for  Optoelectronic  Computing  Systems  has  as  its  specific  goal 
research  on  optical  computing  architectures  suitable  for  use  at  the  highest  possible  speeds. 
The  program  can  be  targeted  toward  exploiting  the  time  domain  because  other  programs 
in  the  Center  are  pursuing  research  on  parallel  optical  systems,  exploiting  optical  intercon¬ 
nection  and  optical  devices  and  materials.  Using  a  general  purpose  computing  architecture 
as  the  focus,  we  are  developing  design  techniques,  tools  and  architectures  for  operation  at 
the  speed  of  light  limit.  Experimental  work  is  being  done  with  the  somewhat  low  speed 
components  currently  available  but  with  architectures  which  will  scale  up  in  speed  as  faster 
devices  developed.  The  design  algorithms  and  tools  developed  for  a  general  purpose,  stored 
program  computer  are  being  applied  to  other  systems  such  as  optically  controlled  optical 
communications  networks. 
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Gropp,  William  D.,  and  David  E.  Keyes:  Domain  decomposition  with  local  mesh  refinement. 
ICASE  Report  No.  91-19,  February  22,  1991,  32  pages.  Submitted  to  SIAM  Journal  of 
Scientific  and  Statistical  Computation. 


We  describe  a  preconditioned  Krylov  iterative  algorithm  based  on  domain  decomposition 
for  linear  systems  arising  from  implicit  finite- difference  or  finite-element  discretizations  of 
partial  differential  equation  problems  requiring  local  mesh  refinement.  To  keep  data  struc¬ 
tures  as  simple  as  possible  for  parallel  computing  applications,  we  define  the  fundamental 
computational  unit  in  the  algorithm  as  a  subregion  of  the  domain  spanned  by  a  locally  uni¬ 
form  tensor-product  grid,  called  a  tile.  In  the  tile-based  domain  decomposition  approach, 
two  levels  of  discretization  are  considered  at  each  point  of  the  domain:  a  global  coarse  grid 
defined  by  tile  vertices  only,  and  a  local  fine  grid  where  the  degree  of  resolution  can  vary  from 
tile  to  tile.  One  global  level  and  one  local  level  provide  the  flexibility  required  to  adaptively 
discretize  a  diverse  collection  of  problems  on  irregular  regions  and  solve  them  at  convergence 
rates  that  deteriorate  only  logarithmically  in  the  finest  mesh  parameter,  with  the  coarse 
tessellation  held  fixed.  A  logarithmic  departure  from  optimality  seems  to  be  a  reasonable 
compromise  for  the  simplicity  of  the  composite  grid  data  structure  and  concomitant  regular 
data  exchange  patterns  in  a  multiprocessor  environment.  We  report  some  experiments  with 
up  to  1024  tiles,  comment  on  the  evolution  of  the  algorithm,  and  contrast  it  with  optimal 
nonrefining  two-level  algorithms  and  optimal  refining  multilevel  algorithms.  Computational 
comparisons  with  some  other  popular  methods  are  presented. 


Gropp,  William  D.,  and  David  E.  Keyes:  Domain  decomposition  methods  in  computational 
fluid  dynamics.  ICASE  Report  No.  91-20,  February  22,  1991,  24  pages.  Submitted  to  Inter¬ 
national  Journal  of  Numerical  Methods  in  Fluids. 


The  divide-and-conquer  paradigm  of  iterative  domain  decomposition,  or  substructuring, 
has  become  a  practical  tool  in  computational  fluid  dynamics  applications  because  of  its 
flexibility  in  accommodating  adaptive  refinement  through  locally  uniform  (or  quasi-uniform) 
grids,  its  ability  to  exploit  multiple  discretizations  of  the  operator  equations,  and  the  modular 
pathway  it  provides  towards  parallelism.  We  illustrate  these  features  on  the  classic  model 
problem  of  flow  over  a  backstep  using  Newton’s  method  as  the  nonlinear  iteration.  Multiple 
discretizations  (second-order  in  the  operator  and  first-order  in  the  preconditioner)  and  locally 
uniform  mesh  refinement  pay  dividends  separately,  and  they  can  be  combined  synergistically. 
We  include  sample  performance  results  from  an  Intel  iPSC/860  hypercube  implementation. 


Tanveer,  S.:  Singularities  in  water  waves  and  Rayleigh- Taylor  instability.  ICASE  Report  No. 
91-21.  February  15,  1991,  39  pages.  Submitted  to  Proc.  Roy.  Soc.  LONDON  A. 


This  paper  is  concerned  with  singularities  in  inviscid  two  dimensional  finite  amplitude 
water  waves  and  inviscid  Rayleigh- Taylor  instability.  For  the  deep  water  gravity  waves  of 
permanent  form,  through  a  combination  of  analytical  and  numerical  methods,  we  present 
results  describing  the  precise  form,  number  and  location  of  singularities  in  the  unphysical 
domain  as  the  wave  height  is  increased.  We  then  show  how  the  information  on  the  singularity 
can  be  used  to  calculate  water  waves  numerically  in  a  relatively  efficient  fashion.  We  also 
show  that  for  two  dimensional  water  waves  in  a  finite  depth  channel,  the  nearest  singularity 
in  the  unphysical  region  has  the  same  form  as  for  deep  water  waves.  However,  associated 
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with  such  a  singularity,  there  is  a  series  of  image  singularities  at  increasing  distances  from  the 
physical  plane  with  possibly  different  behavior.  Further,  for  the  Rayleigh- Taylor  problem  of 
motion  of  fluid  over  vacuum,  and  for  the  unsteady  water  wave  problem,  we  derive  integro- 
differential  equations  valid  in  the  unphysical  region  and  show  how  these  equations  can  give 
information  on  the  nature  of  singularities  for  arbitrary  initial  condiiions. 


Smyrlis,  Yiorgos  S.,  and  Demetrious  T.  Papageorgiou:  Predicting  chaos  for  infinite  dimen¬ 
sional  dynamical  systems:  The  Kuramoto-Sivashinsky  equation,  a  case  study.  ICASE  Report 
No.  91-22,  February  19, 1991,  14  pages.  Submitted  to  Proceedings  of  the  National  Academy 
of  Sciences,  USA. 


The  results  of  extensive  computations  are  presented  in  order  to  accurately  charaterize 
transitions  to  chaos  for  the  Kuramoto-Sivashinsky  equation.  In  particular  we  follow  the 
oscillatory  dynamics  in  a  window  that  supports  a  complete  sequence  of  period  doubling 
bifurcations  preceding  chaos.  As  many  as  thirteen  period  doublings  are  followed  and  used 
to  compute  the  Feigenbaum  number  for  the  cascade  and  so  enable,  for  the  first  time,  an 
accurate  numerical  evaluation  of  the  theory  of  universal  behavior  of  nonlinear  systems,  for 
an  infinite  dimensional  dynamical  system.  Further  more,  the  dynamics  at  the  threshold  of 
chaos  exhibit  a  fractal  behavior  which  is  demonstrated  and  used  to  compute  a  universal 
scaling  factor  that  enables  the  self-similar  continuation  of  the  solution  into  a  chaotic  regime. 


Thangam,  S.,  and  C.G.  Speziale:  Turbulent  separated  flow  past  a  backward-facing  step:  A 
critical  evaluation  of  two-equation  turbulence  models .  ICASE  Report  No.  91-23,  February 
20,  1991,  27  pages.  Submitted  to  Physics  of  Fluids  A. 


The  ability  of  two-equation  turbulence  models  to  accurately  predict  separated  flows  is 
analyzed  from  a  combined  theoretical  and  computational  standpoint.  Turbulent  flow  past  a 
backward  facing  step  is  chosen  as  a  test  case  in  an  effort  to  resolve  the  variety  of  conflicting 
results  that  have  been  published  during  the  past  decade  concerning  the  performance  of  two- 
equation  models.  It  is  found  that  the  errors  in  the  reported  predictions  of  the  K  —  e  model 
have  two  major  origins:  (1)  numerical  problems  arising  from  inadequate  resolution,  and 
(2)  inaccurate  predictions  for  normal  Reynolds  stress  differences  arising  from  the  use  of  an 
isotropic  eddy  viscosity.  Inadequacies  in  near  wall  modeling  play  a  substantially  smaller  role. 
Detailed  calculations  are  presented  which  strongly  indicate  that  the  standard  K  —  e  -  when 
modified  with  an  independently  calibrated  anisotropic  eddy  viscosity  -  can  yield  surprisingly 
good  predictions  for  the  backstep  problem. 


Jacobs,  P.A.,  and  R.J.  Stalker:  Mach  4  and  Mach  8  axisymmetric  nozzles  for  a  shock  tunnel. 
ICASE  Report  No.  91-24,  February  26, 1991,  26  pages.  To  be  submitted  to  the  Aeronautical 
Journal. 


This  study  examines  the  performance  of  two  axisymmetric  nozzles  which  were  designed 
to  produce  uniform,  parallel  flow  with  nominal  Mach  numbers  of  4  and  8.  A  free-piston- 
driven  shock  tube  was  used  to  supply  the  nozzle  with  high-temperature,  high-pressure  test 
gas.  The  inviscid  design  procedure  treated  the  nozzle  expansion  in  two  stages.  Close  to 
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the  nozzle  throat,  the  nozzle  wall  was  specified  as  conical  and  the  gas  flow  was  treated  as 
a  quasi-one-dimensional  chemically-reacting  flow.  At  the  end  of  the  conical  expansion,  the 
gas  was  assumed  to  be  calorically  perfect  and  a  contoured  wall  was  designed  (using  Method- 
of-Characteristics)  to  convert  the  source  flow  into  a  uniform  and  parallel  flow  at  the  end  of 
the  nozzle.  Performance  was  assessed  by  measuring  Pitot  pressures  across  the  exit  plane 
of  the  nozzles  and,  over  the  range  of  operating  conditions  examined,  the  nozzles  produced 
satisfactory  test  flows.  However,  there  were  flow  disturbances  in  the  Mach  8  nozzle  flow  that 
persisted  for  significant  times  after  flow  initiation. 


Mavriplis,  Dimitri  J.:  Unstructured  and  adaptive  mesh  generation  for  high  Reynolds  number 
viscous  flows.  ICASE  Report  No.  91-25,  February  27, 1991,  25  pages.  Proceedings  of  the  3rd 
International  Conference  on  Numerical  Grid  Generation  Conference  to  be  held  in  Barcelona, 
SPAIN,  June  3-7,  1991. 


A  method  for  generating  and  adaptively  refining  a  highly  stretched  unstructured  mesh, 
suitable  for  the  computation  of  high- Reynolds- number  viscous  flows  about  arbitrary  two- 
dimensional  geometries  has  been  developed.  The  method  is  based  on  the  Delaunay  triangu¬ 
lation  of  a  predetermined  set  of  points  and  employs  a  local  mapping  in  order  to  achieve  the 
high  stretching  rates  required  in  the  boundary-layer  wake  regions.  The  initial  mesh-point 
distribution  is  determined  in  a  geometry- adaptive  manner  which  clusters  points  in  regions 
of  high  curvature  and  sharp  corners.  Adaptive  mesh  refinement  is  achieved  by  adding  new 
points  in  regions  of  large  flow  gradients,  and  locally  retriangulating,  thus  obviating  the  need 
for  global  mesh  regeneration.  Initial  and  adapted  meshes  about  complex  multi-element  airfoil 
geometries  are  shown  and  compressible  flow  solutions  are  computed  on  these  meshes. 


Cai,  Wei,  and  Chi- Wang  Shu:  Uniform  high  order  spectral  methods  for  one  and  two  dimen¬ 
sional  Euler  equations.  ICASE  Report  No.  91-26,  March  5,  1991,  52  pages.  Submitted  to 
Journal  of  Computational  Physics. 


In  this  paper  we  study  uniform  high  order  spectral  methods  to  solve  multi-dimensional 
Euler  equations  for  gas  dynamics.  Uniform  high  order  spectral  approximations  with  spectral 
accuracy  in  smooth  regions  of  solutions  are  constructed  by  introducing  the  idea  of  the  Essen¬ 
tially  Non-Oscillatory  (ENO)  polynomial  interpolations  into  the  spectral  methods.  Based 
on  the  new  approximations,  we  propose  nonoscillatory  spectral  methods  which  possess  the 
properties  of  both  upwinding  difference  schemes  and  spectral  methods.  We  present  numerical 
results  for  the  inviscid  Burgers’  equation,  and  for  one  dimensional  Euler  equations  including 
the  interactions  between  a  shock  wave  and  density  disturbance,  Sod’s  and  Lax’s  shock  tube 
problems,  and  the  blast  wave  problem.  Finally,  we  simulate  the  interaction  between  a  Mach 
3  two  dimensional  shock  wave  and  a  rotating  vortex. 


Banks,  H.  T.,  G.  Propst,  and  R.  J.  Silcox:  A  comparison  of  time  domain  boundary  conditions 
for  acoustic  waves  in  wave  guides.  ICASE  Report  No.  91-27. 


We  consider  several  types  of  boundary  conditions  in  the  context  of  time  domain  models 
for  acoustic  waves.  Experiments  with  different  duct  terminations  (hardwall,  free  radiation, 
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foam,  wedge)  were  carried  out  in  a  wave  duct  from  which  reflection  coefficients  over  a  wide 
frequency  range  were  obtained.  These  reflection  coefficients  are  used  to  estimate  parameters 
in  the  time  domain  boundary  conditions  and  a  comparison  of  the  relative  merits  of  the 
models  in  describing  the  data  is  presented. 


Nicolaides,  R.A.:  Analysis  and  convergence  of  the  MAC  Scheme.  1.  The  linear  problem. 
ICASE  Report  No.  91-28,  March  11,  1991,  18  pages.  Submitted  to  SIAM  Journal  of  Nu¬ 
merical  Analysis. 


The  MAC  discretization  of  fluid  flow  is  analyzed  for  the  stationary  Stokes  equations.  It 
is  proved  that  the  discrete  approximations  do  in  fact  converge  to  the  exact  solutions  of  the 
flow  equations.  Estimates  using  mesh  dependent  norms  analogous  to  the  standard  H1  and 
L 2  norms  are  given  for  the  velocity  and  pressure  respectively. 


Sarkar,  S.,  G.  Erlebacher,  and  M.  Y.  Hussaini:  Direct  simulation  of  compressible  turbulence  in 
a  shear  flow.  ICASE  Report  No.  91-29,  March  13, 1991,  41  pages.  Submitted  to  Theoretical 
and  Computational  Fluid  Dynamics. 


The  purpose  of  this  study  is  to  investigate  compressibility  effects  on  the  turbulence  in  ho¬ 
mogeneous  shear  flow.  We  find  that  the  growth  of  the  turbulent  kinetic  energy  decreases  with 
increasing  Mach  number  -  a  phenomenon  which  is  similar  to  the  reduction  of  turbulent  ve¬ 
locity  intensities  observed  in  experiments  on  supersonic  free  shear  layers.  An  examination  of 
the  turbulent  energy  budget  shows  that  both  the  compressible  dissipation  and  the  pressure- 
dilatation  contribute  to  the  decrease  in  the  growth  of  kinetic  energy.  The  pressure-dilatation 
is  predominantly  negative  in  homogeneous  shear  flow,  in  contrast  to  its  predominantly  pos¬ 
itive  behavior  in  isotropic  turbulence.  The  different  signs  of  the  pressure-dilatation  are 
explained  by  theoretical  consideration  of  the  equations  for  the  pressure  variance  and  density 
variance.  We  obtained  previously  the  following  results  for  isotropic  turbulence;  first,  the 
normalized  compressible  dissipation  is  of  O(Mf),  and  second,  there  is  approximate  equipar- 
tition  between  the  kinetic  and  potential  energies  associated  with  the  fluctuating  compressible 
mode.  Both  these  results  have  now  been  substantiated  in  the  case  of  homogeneous  shear. 
The  dilatation  field  is  significantly  more  skewed  and  intermittent  than  the  vorticity  field. 
Strong  compressions  seem  to  be  more  likely  than  strong  expansions. 


Pruett,  C.  David,  Lian  L.  Ng,  and  Gordon  Erlebacher:  On  the  nonlinear  stability  of  a  high¬ 
speed,  axisymmetric  boundary  layer.  ICASE  Report  No.  91-30,  March  19,  1991,  46  pages. 
To  be  submitted  to  Physics  of  Fluids. 


The.  stability  of  a  high-speed,  axisymmetric  boundary  layer  is  investigated  using  sec¬ 
ondary  instability  theory  and  direct  numerical  simulation.  Parametric  studies  based  on  tem¬ 
poral  secondary  instability  theory  indentify  subharmonic  secondary  instability  as  a  likely 
path  to  transition  on  a  cylinder  at  Mach  4.5.  The  theoretical  predictions  are  validated  by 
direct  numerical  simulation  of  temporally-evolving  primary  and  secondary  disturbances  in  an 
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axisymmetric  boundary-layer  flow.  At  small  amplitudes  of  the  secondary  disturbance,  pre¬ 
dicted  growth  rates  agree  to  several  significant  digits  with  values  obtained  from  the  spectrally- 
accurate  solution  of  the  compressible  Navier-Stokes  equations.  Qualitative  agreement  per¬ 
sists  to  large  amplitudes  of  the  secondary  disturbance.  Moderate  transverse  curvature  is 
shown  to  significantly  affect  the  growth  rate  of  axisymmetric  “second  mode”  disturbances, 
the  likely  candidates  of  primary  instability.  The  influence  of  curvature  on  secondary  insta¬ 
bility  is  largely  indirect  but  most  probably  significant,  through  modulation  of  the  primary 
disturbance  amplitude.  Subharmonic  secondary  instability  is  shown  to  be  predominantly 
inviscid  in  nature,  and  to  account  for  spikes  in  the  Reynolds  stress  components  at  or  near 
the  critical  layer. 


Cai,  Wei,  David  Gottlieb,  and  Chi-Wang  Shu:  On  one-sided  filters  for  spectral  Fourier 
approximations  of  discontinuous  functions.  ICASE  Report  No.  91-31,  March  28,  1991,  19 
pages.  Submitted  to  SIAM  Journal  of  Numerical  Analysis. 


In  this  paper  we  prove  the  existence  of  one-sided  filters,  for  spectral  Fourier  approxima¬ 
tions  of  discontinuous  functions,  which  can  recover  spectral  accuracy  up  to  the  discontinuity 
from  one  side.  We  also  use  a  least  square  procedure  to  construct  such  a  filter  and  test  it  on 
several  discontinuous  functions  numerically. 


Cockburn,  Bernardo,  and  Chi-Wang  Shu:  The  P1  —  RKDG  method  for  two-dimensional 
Euler  equations  of  gas  dynamics.  ICASE  Report  No.  91-32,  March  28,  1991,  12  pages.  Sub¬ 
mitted  to  Second  International  Symposium  on  High  Performance  Computing. 


We  continue  our  earlier  work  on  a  class  on  nonlinearly  stable  Runge-Kutta  local  projec¬ 
tion  discontinuous  Galerkin  (RKDG)  finite  element  methods  for  conservation  laws.  Two- 
dimensional  Euler  equations  for  gas  dynamics  are  solved  using  P1  elements.  We  discuss 
the  generalization  of  the  local  projection,  which  for  scalar  nonlinear  conservation  laws  was 
designed  to  satisfy  a  local  maximum  principle,  to  systems  of  conservation  laws  such  as  the 
Euler  equations  of  gas  dynamics  using  local  characteristic  decompositions.  Numerical  exam¬ 
ples  include  the  standard  regular  shock  reflection  problem,  the  forward  facing  step  problem 
and  the  double  Mach  reflection  problem.  These  preliminary  numerical  examples  are  chosen 
to  show  the  capacity  of  our  approach  to  obtain  nonlinearly  stable  results  comparable  with 
the  modern  nonoscillatory  finite  difference  methods.  Generalizations  to  Pfc  elements  with 
k  >  1  and  the  use  of  adaptive  triangulations  to  minimize  local  errors  constitute  ongoing 
research. 


Bryan,  Kurt:  Numerical  recovery  of  certain  discontinuous  electrical  conductivities.  ICASE 
Report  No.  91-33,  March  28,  1991,  28  pages.  Submitted  to  Inverse  Problems. 

The  inverse  problem  of  recovering  an  electrical  conductivity  of  the  form  'y(x)  =  1  +  (k  — 
l)Xi>  (Xd  is  the  characteristic  function  of  D)  on  a  region  Q  C  1R2  from  boundary  data  is 
considered,  where  D  CC  and  k  is  some  positive  constant.  A  linearization  of  the  forward 
problem  is  formed  and  used  in  a  least  squares  output  method  for  approximately  solving  the 
inverse  problem.  Convergence  results  are  proved  and  some  numerical  results  presented. 
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ICASE  INTERIM  REPORTS 


Eidson,  T.  M.,  and  T.  A.  Zang:  Filtering  analysis  of  a  direct  numerical  simulation  of  the 
turbulent  Rayleigh-Benard  problem.  Interim  Report  No.  15,  December  3,  1990,  49  pages. 


A  filtering  analysis  of  a  turbulent  flow  has  been  developed  which  provides  details  of 
the  path  of  tne  kinetic  energy  of  the  flow  from  its  creation  via  thermal  production  to  its 
dissipation.  A  low-pass  spatial  filter  is  used  to  split  the  velocity  and  the  temperature  field 
into  a  filtered  component  (composed  mainly  of  scales  larger  than  a  specific  size,  nominally 
the  filter  width)  and  a  fluctuation  component  (scales  smaller  than  a  specific  size).  Variables 
derived  from  these  fields  can  fall  mto  one  of  the  above  two  ranges  or  be  composed  of  a  mixture 
of  scales  dominated  by  scales  near  the  specific  size.  The  filter  is  used  to  split  the  kinetic 
energy  equation  into  three  equations  corresponding  to  the  three  scale  ranges  described  above. 

The  data  from  a  direct  simulation  of  the  Rayleigh-Benard  problem  for  conditions  where 
the  flow  is  turbulent  is  used  to  calculate  the  individual  terms  in  the  three  kinetic  energy 
equations.  This  is  done  for  a  range  of  filter  widths.  These  results  are  used  to  study  the 
spatial  location  and  the  scale  range  of  the  thermal  energy  production,  the  cascading  of 
kinetic  energy,  the  diffusion  of  kinetic  energy  and  the  energy  dissipation.  These  results 
are  used  also  to  evaluate  two  subgrid  models  typically  used  in  large-eddy  simulations  of 
turbulence.  Subgrid  models  attempt  to  model  the  energy  below  the  filter  width  that  is 
removed  by  a  low-pass  filter. 


Jacobs,  P.A.,  R.G.  Morgan,  R.C.  Rogers,  M.  Wendt,  C.  Brescianini,  A.  Pauli,  and  G.  Kelly: 
Preliminary  calibration  -/  «  generic  scramjet  combustor.  Interim  Report  No,  16,  March  22, 
1991,  31  pages. 


The  results  of  a  preliminary  investigation  of  the  combustion  of  hydrogen  fuel  at  hypersonic 
flow  conditions  are  provided.  The  tests  were  performed  in  a  generic,  constant-area  combustor 
model  with  test  gas  supplied  by  a  free-piston-driven  reflected-shock  tunnel.  Static  pressure 
measurements  along  the  combustor  wall  indicated  that  burning  did  occur  for  combustor  inlet 
conditions  of  PataUc  -  19 kPa,  Tatatic  —  1080 A  and  U  ~  3630m/s  with  a  fuel  equivalence 
ratio  (j>  ~  0.9.  These  inlet  conditions  were  obtained  by  operating  the  tunnel  with  stagnation 
enthalpy  H,  8.1  MJ/kg,  stagnation  pressure  P„  ~  52 MPa  and  a  contoured  nozzle  with  a 
nominal  exit  Mach  number  of  5.5. 
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ICASE  COLLOQUIA 

October  1,  1990  through  March  31,  1991 


Name  /Affiliation /Title 


Date 


Dr.  Saul  Abarbanel,  Tel-Aviv  University  and  ICASE  October  3 

“Spurious  Frequencies  as  a  Result 
of  Numerical  Boundary  Treatments” 

Dr.  Rahul  Simha,  The  College  of  William  and  Mary  October  12 

“Optimization  of  Resource  Control  in  Commmunication  Systems” 

M.  Zubair,  Old  Dominion  University  October  16 

“Parallel  Partitioning  of  Sparse  Matrices” 

Patrick  deBondeli,  CR2A,  France 

and  Dennis  B.  Mulcare,  LASC-Georgia  October  16 

“Concurrent  Embedded  System  Architecture  Design  Methods” 

Dr.  Eli  Reshotko,  Case  Western  Reserve  University  October  18 

“Effect  of  Elevated  Free- Stream  Turbulence 
on  Turbulent  Skin  Friction” 

Vidyadhar  Y.  Mudkavi,  California  Institute  of  Technology  October  22 

“Axisymmetric  Waves  in  Vortex  Filaments” 

Dr.  Richard  Barnwell,  NASA  Langley  Research  Center  October  23 

“Nonadiabatic  and  Three-Dimensional  Effects  in  Compressible 
Turbulent  Boundary  Layers” 

Cathy  Westbury,  Queen’s  University,  Ontario,  Canada  October  24 

“Turbulent  Cooling  Flows  in  Clusters  of  Galaxies” 

Professor  Harry  Jordan,  University  of  Colorado  November  1 

“Time-Space  Trade  Offs  in  Optical  Computing” 

Andy  Kopser,  Supercomputing  Research  Center  November  2 

“SPLASH:  A  Linear  Reconfigurable  Logic  Array” 

Dr.  Randolph  A.  Graves,  Supercomputing  Systems,  Inc.  November  5 

“Distributed  Memory  Parallel  Processing 
with  High-Performance  Communications” 
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Name/Affiliation/Title 


Date 


Dr.  Sedat  Biringen,  University  of  Colorado 
“Numerical  Simulations  of  Incompressible 
and  Compressible  Flow  Stability” 

Professor  Lotfi  Zadeh,  University  of  California,  Berkeley 
“Fuzzy  Logic  and  its  Applications” 

Dr.  Yuh-Roung  Ou,  ICASE 
“Control  of  a  Cylinder  Wake  by  Rotation” 

Dr.  Linda  Krai,  McDonnell  Douglas  Research  Laboratories 
“Computational  Studies  in  Transition  and  Turbulence” 

Dr.  Hassan  Aref,  University  of  California  at  San  Diego 
“Lagrangian  CFD” 

Steven  Ashby,  Lawrence  Livermore  National  Laboratory 
“A  Matrix  Analysis  of  Conjugate  Gradient  Algorithms” 

Dr.  G.  Comte-Bellot,  Ecole  Centrale  de  Lyon,  France 
“Propagation  of  Acoustic  Waves  Through  Turbulence” 

Professor  A.  M.  Yaglom,  Institute  of  Atmospheric  Physics, 

USSR  Academy  of  Sciences 
“Spectra  and  Correlation  Functions  of  Turbulence 
in  any  Unstably  Stratified  Atmospheric  Surface  Layer” 

Professor  H.  Peerhossaini,  University  of  Nantes  -  France 
“Boundary  Layer  Flow  on  the  Turbine  Blade  Pressure  Side 
(Gortler  Vortex  Flow)” 

Dr.  Yiorgos  Smyrlis,  University  of  California,  Los  Angeles 
“The  Route  to  Chaos  for  the  Kuramoto-Sivaskinsky  Equation” 

Dr.  Craig  Douglas,  Yale  University 
“A  Greatly  Simplified  Theory  for  Multigrid  Methods” 

Dr.  Alvin  Bayliss,  Northwestern  University 
“Nonlinear  Dynamics  in  Combustion” 

Donald  Delisi,  Northwest  Research  Associated,  Inc. 

“New  Laboratory  Measurements  to  the  Evolution  of  a  Vortex 
Pair:  Is  CFD  up  to  the  Challenge?” 


November  5 


November  6 


November  6 


November  ? 


November  9 


November  13 


November  14 


November  14 


November  15 


November  21 


November  30 


December  4 


December  5 
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Name/Affiliation/Title 


Date 


Professor  Ralph  Metcalfe,  University  of  Houston 
“Large  Scale  Structures  in  Reacting  Mixing  Layers” 

Dr.  David  Kopriva,  The  Florida  State  University 
“Spectral  Solution  of  Supersonic  Flows  over  Cones  and  Wedges” 

Dr.  Phillip  Hall,  University  of  Manchester,  England 
“Vortex  Instabilities  in  Hypersonic  Flows: 

Sutherland  Law  Fluids  and  Real  Gas  Effects” 

Professor  Z.U.A.  Warsi,  Mississippi  State  University 
“Numerical  Mapping  of  Arbitrary  Domains  Using  Spectral  Methods” 

Professor  Anthony  P.  Reeves,  Cornell  University 
“Programming  Highly  Parallel  Multicomputers: 

The  Paragon  Approach” 

Xinming  A.  Lin,  University  of  California,  Los  Angeles 
“An  Object-Oriented  Particle  Simulation  Approach” 

Dr.  Siva  M.  Mangalam,  Analytical  Services  &  Materials,  Inc. 

“Some  Advances  in  Spatial  and  Temporal  Flow  Diagnostics 
Techniques  for  Flight,  Wind  Tunnel,  and  CFC  Research” 

Jurgen  Schoene,  German  Aerospace  Research  Establishment 
“Design  of  Supersonic  Wings  Using  an  Optimization  Strategy 
Coupled  with  a  Solution  Scheme  for  the  Euler  Equations” 

Dr.  N.  Knoll,  German  Aerospace  Research  Establishment 
“A  Quadrilateral  Based  Second-Order  TVD  Method 
for  Unstructured  Adaptive  Meshes” 

Dr,  Reda  R.  Mankbadi,  NASA  Lewis  Research  Center 
“Boundary-Layer  Transition:  Developing  a  Critical-Layer  Theory” 

Dr.  Ugo  Piomelli,  University  of  Maryland,  College  Park 
“A  Dynamic  Subgrid  Scale  Stress  Model  for  Large- Eddy  Simulations” 

Roldan  Pozo,  University  of  Colorado  at  Boulder 
“Performance  Analysis  of  Parallel  Architectures 
for  Scientific  Computing” 

John  Van  Rosendale,  ICASE 
“An  Optimal  Parallel  Time  Multigrid  Algorithm” 
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December  6 

December  10 

December  12 

January  4 
January  9 

January  11 

January  16 

January  18 

January  18 

January  23 
February  15 
February  20 

February  22 


Name  /  Affiliation /Title 


Date 


Ammar  Degani,  Lehigh  University 
'‘The  Asymptotic  Structure  of  a  Three-Dimensional  Turbulent 
Boundary  Layer  -  Theory  and  Application” 

Professor  Vladimar  Rokhlin,  Yale  University 
“Fast  Evaluation  of  Functions  of  Dense  Matrices” 

Rupak  Biswas,  Renssalaer  Polytechnic  Institute 
“Parallel  Adaptive  Methods  for  Hyperbolic  Partial 
Differential  Equations” 

Phil  Keenan,  University  of  Chicago 
“Thermal  Simulation  of  Pipeline  Flow” 

James  Quirk,  Cranfield  Institute  of  Technology 
“An  Adaptive  Grid  Algorithm  for  Computational  Shock  Hydrodynamics” 

Dinshaw  Balsara,  The  Johns  Hopkins  University 
“Jet  Disruption  Using  High-Resolution  Godunov  Schemes” 

Professor  Alok  N.  Choudhury,  Syracuse  University 
“FORTRAN77D  and  FORTRAN90D:  Scalable  and  Portable  Software 
for  SIMD  and  MIMD  Distributed  Memory  Parallel  Computers” 

David  Kotz,  Duke  University 
“Prefetching  and  Caching  Techniques  in  File  Systems 
for  MIMD  Multiprocessors” 

Professor  David  G.  Crighton,  University  of  Cambridge,  England 
“The  Noise  of  a  Modern  Many-Bladed  Propeller-Predictions 
Based  on  Asymptotic  Theory” 

Paul  Durbin,  Stanford/ Ames  Center  for  Turbulence  Research 
“Turbulence  Modeling  Near  Solid  Boundaries” 


February  28 

March  1 

March  7 

March  8 

March  11 

March  12 

March  14 

March  18 

March  20 

March  28 
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OTHER  ACTIVITIES 


A  workshop  on  “Unstructured  Scientific  Computation  on  Scalable  Multiprocessors”  spon¬ 
sored  by  ICASE  was  held  at  the  Ramada  Inn,  Kill  Devil  Hills,  NC  on  October  29-31,  1990. 
Sixty-five  people  attended.  The  objective  of  this  workshop  was  to  discuss  mathematical  and 
computer  science  issues  related  to  unstructured  scientific  computation  on  multiprocessor 
systems.  A  proceedings  will  be  published  by  the  M.I.T.  Press. 

On  December  17-19  ICASE  hosted  a  workshop  to  identify  research  opportunities  in  soft¬ 
ware  of  importance  to  NASA’s  initiative  under  the  Federal  High  Performance  Computing 
and  Communications  (HPCC)  program.  Approximately  40  people  participated  representing 
the  NASA  centers  involved  in  HPCC  and  several  universities.  As  a  result  of  the  workshop 
a  matrix  was  created  that  linked  important  areas  of  research  in  software  with  organizations 
that  were  1  epared  to  pursue  the  research. 

On  November  12-13,  1990  a  symposium  on  ’’Recent  Development  in  Turbulence”  was 
sponsored  by  ICASE  and  NASA  Langley  Research  Center  was  held  at  the  Omni  Hotel, 
Newport  News,  VA.  Forty-eight  people  attended  this  syir*  osium.  The  objectives  were  to 
review  the  latest  thrusts  in  turbulence  research,  including  experimental  and  theoretical  re¬ 
search.  Recent  work  in  turbulence  modeling,  as  well  as  current  work  on  the  development 
of  the  Proper  Orthogonal  Decomposition  method,  was  also  discussed.  Papers  from  this 
symposium  will  be  published  by  Springer- Verlag  in  a  volume  entitled  Studies  in  Turbulence. 
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ICASE  STAFF 


I.  ADMINISTRATIVE 

Robert  G.  Voigt,  Director  Ph.D.,  Mathematics,  University  of  Maryland,  1969  Numerical  and 
Algorithms  for  Parallel  Computers 

Linda  T.  Johnson,  Office  and  Financial  Administrator 

Etta  M.  Blair,  Personnel/ Bookkeeping  Secretary 

Barbara  A.  Cardasis,  Administrative  Secretary 

Rosa  H.  Milby,  Technical  Publications/Summer  Housing  Secretary 

Shelly  D.  Millen,  Technical  Publications  Secretary 

Emily  N.  Todd,  Executive  Secretary/Visitor  Coordinator 

II.  SCIENCE  COUNCIL  for  APPLIED  MATHEMATICS  and  COMPUTER 
SCIENCE 

Tony  Chen,  Professor,  Department  of  Mathematics,  University  of  California  at  Los  Angeles. 
S.  Lennart  Johnsson,  Thinking  Machines  Corporation. 

Michael  O’Donnell,  Professor,  Department  of  Computer  Science,  University  of  Chicago. 

Robert  O’Malley,  Jr.,  Chairman,  Department  of  Mathematical  Sciences,  Rensselaer  Poly¬ 
technic  Institute. 

Joseph  Oliger,  Professor,  Computer  Science  Department,  Stanford  University. 

Stanley  J.  Osher,  Professor,  Mathematics  Department,  University  of  California. 

John  Rice,  Chairman,  Department  of  Computer  Science,  Purdue  University. 

Burton  Smith,  Tera  Computer  Company,  Seattle,  WA. 

Robert  G.  Voigt,  Director,  Institute  for  Computer  Applications  in  Science  and  Engineering, 
NASA  Langley  Research  Center. 

HI.  ASSOCIATE  MEMBERS 

Saul  S.  Abarbanel,  Professor,  Department  of  Applied  Mathematics,  Tel-Aviv  University, 
Israel. 
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H.  Thomas  Banks,  Professor,  Center  for  Applied  Mathematical  Sciences,  University  of  South¬ 
ern  California. 

David  Gottlieb,  Professor,  Division  of  Applied  Mathematics,  Brown  University. 

Peter  D.  Lax,  Professor,  Courant  Institute  of  Mathematical  Sciences,  New  York  University. 
Merrell  L.  Patrick,  New  Technologies  Program  Director,  National  Science  Foundation. 


IV.  CHIEF  SCIENTIST 

M.  Yousuff  Hussaini  -  Ph.D.,  Mechanical  Engineering,  University  of  California,  1970.  Com¬ 
putational  Fluid  Dynamics.  (Beginning  April  1978) 


V.  LEAD  COMPUTER  SCIENTIST 

Joel  H.  Saltz  -  Ph.D.,  Computer  Science,  Duke  University,  1985.  Parallel  Computing  with 
Emphasis  on  Systems  and  Algorithmic  Issues.  (Beginning  July  1989) 


V.  SENIOR  STAFF  SCIENTIST 

Gordon  Erlebacher  -  Ph.D.,  Plasma  Physics,  Columbia  University,  1983.  Computational 
Fluid  Dynamics.  (November  1989  to  November  1994) 

Dimitri  Mavriplis  -  Ph.D.,  Mechanical  and  Aerospace  Engineering,  Princeton  University, 
1988  Grid  Techniques  for  Computational  Fluid  Dynamics.  (February  1987  to  September 
1985) 

Charles  G.  Speziale  -  Ph.D.,  Aerospace  and  Mechanical  Sciences,  Princeton  University, 
1978.  Fluid  Dynamics  with  Emphasis  on  Turbulence  Modeling  and  the  Transition  Process. 
(September  1987  to  September  1991) 

John  R.  Van  Rosendale  -  Ph.D.,  Computer  Science,  University  of  Illinois,  1980.  Parallel 
Systems  and  Algorithms.  (July  1,  1989  to  July  1,  1992) 


VI.  SCIENTIFIC  STAFF 

Kurt  M.  Bryan  -  Ph.D.,  Mathematics,  University  of  Washington,  1990.  Theoretical  and 
Computational  Methods  for  Inverse  Problems.  (August  1990  to  August  1992) 

Richard  G.  Carter  -  Ph.D.,  Numerical  Analysis,  Rice  University,  1986.  Numerical  Methods 
for  Optimization  Problems.  (September  1987  to  December  1990) 

Leon  M.  Clancy  -  B.S.,  Mechanical  Engineering,  University  of  Washington,  1971.  System 
Manager.  (December  1989  to  Present) 

Thomas  W.  Crockett  -  B.S.,  Mathematics,  College  of  William  and  Mary,  1977.  Parallel 
Systems  Research.  (February  1987  to  February  1992) 

Thomas  M.  Eidson  -  Ph.D.,  Mechanical  Engineering,  University  of  Michigan,  1982.  Parallel 
Techniques  for  Computational  Fluid  Dynamics.  (August  1989  to  August  1991) 
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Yuh-Roung  Ou  -  Ph.D.,  Aerospace  Engineering,  University  of  Southern  California,  1988. 
Control  Systems  for  Fluid  Dynamics.  (November  1988  to  September  1991) 

Peter  W.  Protzel  -  Ph.D.,  Electrical  Engineering,  Technical  University  of  Braunschweig, 
Germany,  1987.  Reliability  of  Computing  Systems.  (March  1987  to  September  1991) 

Sutanu  Sarkar  -  Ph.D.,  Mechanical  and  Aerospace  Engineering,  Cornell  University,  1988. 
Fluid  Dynamics,  Turbulence  Modeling,  Compressible  Turbulence.  (September  1988  to  Septem¬ 
ber  1991) 

Jeffrey  S.  Scroggs  -  Ph.D.,  Computer  Science,  University  of  Illinois  at  Urbana,  1988.  Domain 
Decomposition  for  Differential  Equations.  (July  1988  to  January  1991) 

Sharon  0.  Seddougui  -  Ph.D.,  Applied  Mathematics,  University  of  Exeter,  England,  1988. 
Compressible  Fluid  Dynamics.  (June  1988  to  August  1991) 

Ralph  C.  Smith  -  Ph.D.,  Numerical  Analysis,  Montana  State  University,  1990.  Theoretical 
and  Computational  Issues  Associated  with  Inverse  Problems.  (August  1990  to  August  1992) 


VII.  VISITING  SCIENTISTS 

Shahid  Bokhari  -  Ph.D.,  Electrical  and  Computer  Engineering,  University  of  Massachusetts, 
Amherst,  1978.  Professor,  Electrical  Engineering,  University  of  Engineering  k  Technology, 
Lahore,  Pakistan.  Parallel  Computing  Systems.  (October  1989  -  January  1991) 

Peter  W.  Duck  -  Ph.D.,  Fluid  Mechanics,  University  of  Southampton,  United  Kingdom, 
1975.  Lecturer  in  Mathematics,  Department  of  Mathematics,  University  of  Manchester, 
United  Kingdom.  Numerical  Solution  of  Unsteady  Boundary  Layer  Equations.  (January 
1991) 

Kathleen  E.  Crowley  -  Ph.D.,  Computer  Science,  University  of  Washington,  1986.  Sr.  Pro¬ 
grammer,  Department  of  Computer  Science,  Yale  University.  Environments  for  Parallel 
Computing  Systems.  (February  to  November  1991) 

David  E.  Keyes  -  Ph.D.,  Applied  Mathematics,  Harvard  University,  1984.  Associate  Profes¬ 
sor,  Department  of  Mechanical  Engineering,  Yale  University.  Parallelization  of  Numerical 
Procedures  for  Combustion.  (August  to  December  1990) 

Peter  A.  Jacobs  -  Ph.D.,  Computational  Fluid  Dynamics,  University  of  Queensland,  Aus¬ 
tralia,  1987.  Compressible  Fluid  Dynamics.  (Januray  1990  to  January  1992) 

Piyush  Mehrotra  -  Ph.D.,  Computer  Science,  University  of  Virginia,  1982.  Assistant  Pro¬ 
fessor,  Department  of  Computer  Science,  Purdue  University.  Programming  Languages  for 
Multiprocessor  Systems.  (January  to  September  1991) 

Eitan  Tadmor  -  Ph.D.,  Numerical  Analysis,  Tel-Aviv  University,  1979.  Professor,  Depart¬ 
ment  of  Applied  Mathematics,  Tel-Aviv  University.  Numerical  Methods  for  Partial  Differ¬ 
ential  Equations.  (February  1991) 

Siva  Thangam  -  Ph.D.,  Mechanical  Engineering,  Rutgers  University,  1980.  Professor,  De¬ 
partment  of  Mechanical  Engineering,  Stevens  Institute  of  Technology.  Computational  Fluid 
Mechanics.  (August  1990  to  August  1991) 
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VIII.  CONSULTANTS 


Loyce  M.  Adams  -  Ph.D.,  Applied  Mathematics,  University  of  Virginia,  1983.  Associate  Pro¬ 
fessor,  Department  of  Applied  Mathematics,  University  of  Washington.  Numerical  Methods 
for  Parallel  Computing  Systems. 

Alvin  Bayliss  -  Ph.D.,  Mathematics,  New  York  University,  1975.  Associate  Professor,  Tech¬ 
nological  Institute,  Northwestern  University.  Numerical  Methods  for  Partial  Differential 
Equations. 

Dinshaw  Balsara  -  Ph.D.,  Computational  Fluid  Dynamics,  Astro  Physics,  University  of 
Illinois  at  Urbana,  1990.  Department  of  Physics  and  Astronomy,  Johns  Hopkins  University. 
Parallel  Implementation  on  Adaptive  Godunov  Schemes. 

Marsha  J.  Berger  -  Ph.D.,  Numerical  Analysis,  Stanford  University,  1982.  Research  Asso¬ 
ciate,  Courant  Institute  of  Mathematical  Sciences.  Numerical  Methods  for  Partial  Differen¬ 
tial  Equations. 

Harry  Scott  Berryman  -  B.S.,  Computer  Science,  Yale  University,  1988.  Graduate  Fellow, 
Department  of  Computer  Science,  Yale  University.  Software  for  Parallel  Systems. 

Percy  Bobbitt  -  B.S.,  Aeronautics,  Catholic  University  of  America,  1949.  NASA  Langley 
Research  Center  -  Retired.  Fluid  Mechanics. 

Achi  Brandt  -  Ph.D.,  Mathematics,  Weizmann  Institute  of  Science,  1965.  Professor,  Applied 
Mathematics  Department,  Weizmann  Institute  of  Science,  Israel.  Multigrid  Methods. 

John  A.  Burns  -  Ph.D.,  Mathematics,  University  of  Oklahoma,  1973.  Professor,  Virginia 
Polytechnic  Institute  and  State  University.  Numerical  Methods  in  Feedback  Control  and 
Parameter  Estimation. 

Peter  R.  Eiseman  -  Ph.D.,  Mathematics,  University  of  Illinois,  1970.  Senior  Research  Sci¬ 
entist  and  Adjunct  Professor,  Department  of  Applied  Physics  and  of  Nuclear  Engineering, 
Columbia  University.  Computational  Fluid  Dynamics. 

Robert  E.  Fennell  -  Ph.D.,  Mathematics,  University  of  Iowa,  1969.  Professor,  Department 
of  Mathematical  Sciences,  Clemson  University.  Control  Theory  for  Multivariable  Systems. 

Joel  H.  Ferziger  -  Ph.D.,  Nuclear  Engineering,  University  of  Michigan,  1962.  Professor, 
Thermosciences  Division,  Department  of  Mechanical  Engineering,  Stanford  University.  Fluid 
Dynamics. 

James  F.  Geer  -  Ph.D.,  Applied  Mathematics,  New  York  University,  1967.  Professor,  Sys¬ 
tems  Science  and  Mathematical  Sciences,  Watson  School  of  Engineering,  Applied  Science 
and  Technology,  SUNY-Binghamton.  Perturbation  Methods  and  Asymptotic  Expansions  of 
Solutions  to  Partial  Differential  Equations. 

Chester  E.  Grosch  -  Ph.D.,  Physics  -  Fluid  Dynamics,  Stevens  Institute  of  Technology, 
1967.  Professor,  Department  of  Computer  Science  and  Slover  Professor,  Department  of 
Oceanography,  Old  Dominion  University.  Hydrodynamic  Stability,  Computational  Fluid 
Dynamics,  Unsteady  Boundary  Layers  and  Algorithms  for  Array  Processors. 

Philip  Hall  -  Ph.D.,  Mathematics,  Imperial  College,  England,  1973.  Professor,  Department 
of  Mathematics,  University  of  Manchester,  England.  Computational  Fluid  Dynamics. 
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Amiram  Harten  -  Ph.D.,  Mathematics,  New  York  University,  1974.  Associate  Professor, 
Department  of  Mathematics,  Tel-Aviv  University,  Israel.  Numerical  Solution  for  Partial 
Differential  Equations. 

Thorwald  Herbert  -  Ph.D.,  Aerospace  Engineering,  University  of  Stuttgart,  Germany  1.978. 
Professor,  Department  of  Mechanical  Engineering,  Ohio  State  University.  Fluid  Dynamics. 

Kazufumi  Ito  -  Ph.D.,  Systems  Science  and  Mathematics,  Washington  University,  1981. 
Assistant  Professor,  Department  of  Mathematics,  University  of  Southern  California.  Control 
Theory. 

Thomas  L.  Jackson  -  Ph.D.,  Mathematics,  Rensselaer  Polytechnic  Institute,  1985.  Assistant 
Professor,  Department  of  Mathematics,  Old  Dominion  University.  Numerical  and  Analytical 
Methods  for  Chemically  Reacting  Flows. 

Antony  Jameson  -  Ph.D.,  Magnetohydro-Dynamics,  Cambridge  University,  England,  1963. 
James  S.  McDonnell  Distinguished  Professor,  Department  of  Mechanical  and  Aerospace 
Engineering,  Princeton  University.  Computational  Fluid  Dynamics. 

Charles  R.  Johnson  -  Ph.D.,  Mathematics  and  Economics,  California  Institute  of  Technology, 
1972.  Professor,  Department  of  Mathematics,  College  of  William  and  Mary.  Numerical 
Linear  Algebra. 

Mark  T.  Jones  -  Ph.D.,  Computer  Science,  Duke  University,  1990.  Assistant  Computer 
Scientist,  MCS  Division,  Argonne  National  Labs.  Parallel  Algorithms  for  Numerical  Linear 
Algebra. 

Harry  F.  Jordan  -  Ph.D.,  Physics,  University  of  Illinois,  1967.  Professor  Department  of  Elec¬ 
trical  and  Computer  Engineering,  University  of  Colorado  at  Boulder.  Parallel  Computation. 

Ashwani  K.  Kapila  -  Ph.D.,  Theoretical  and  Applied  Mechanics,  Cornell  University,  1975. 
Associate  Professor,  Department  of  Mathematical  Sciences,  Rensselaer  Polytechnic  Institute. 
Ordinary  and  Partial  Differential  Equations,  Asymptotic  Methods. 

Edward  J.  Kerschen  -  Ph.D.,  Mechanical  Engineering,  Stanford  University,  1978.  Associate 
Professor,  Department  of  Aerospace  and  Mechanical  Engineering,  University  of  Arizona. 
Flow  Dynamics. 

David  E.  Keyes  -  Ph.D.,  Applied  Mathematics,  Harvard  University,  1984.  Assistant  Pro¬ 
fessor,  Mechanical  Engineering,  Yale  University.  Parallelizaton  of  Numerical  Procedures 
Appropriate  for  the  Study  of  Combustion. 

Fumio  Kojima  -  Ph.D.,  Control  Theory,  Kyoto  University,  Japan,  1985.  Visiting  Research 
Assistant  Professor,  Center  for  Applied  Mathematical  Sciences,  University  of  Southren  Cal¬ 
ifornia.  Probabilistic  and  Stochastic  Methods  for  Optimal  Control  Problems. 

Heinz-Otto  Kreiss  -  Ph.D.,  Mathematics,  Royal  Institute  of  Technology,  Sweden  1960.  Pro¬ 
fessor,  Department  of  Applied  Mathematics,  California  Institute  of  Technology.  Numerical 
Analysis. 

William  D.  Lakin  -  Ph.D.,  Applied  Mathematics,  University  of  Chicago,  1968.  Professor, 
Department  of  Mathematics  and  Statistics,  University  of  Vermont.  Computational  Fluid 
Dynamics. 
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